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I.  INTRODUCTION 


This  thesis  has  two  objectives.  The  first  is  to  resolve  discrepan¬ 
cies  in  previously  measured  cross  sections  for  the  ^B(p,n)^C  reaction. 
The  second  is  to  determine  the  extent  to  which  this  reaction  can  be 
used  to  measure  boron  concentrations  in  solids. 

The  ^B(p,n)^C  reaction  was  studied  using  protons  in  the  energy 
range  from  threshold  (3.016  MeV)  to  4.6  MeV.  In  this  range  only  the 
ground  state  of  can  be  produced  by  this  reaction.  The  ^B(p,n)^C 

reaction  has  a  nuclear  disintegration  energy  (Q  value)  of  -2.763  MeV 

12  1 

and  forms  a  compound  nucleus  of  C.  A  partial  energy  diagram  of  the 

12 

compound  nucleus  is  shown  in  Figure  1.  Two  other  decay  modes  of  C, 

3 

which  are  not  studied  here,  are:  the  production  of  Be  +  a,  with  a 

12 

center  of  mass  energy  of  7.367  MeV  with  respect  to  the  ground  state  C 

12 

compound  nucleus,  and  C  +  y. 

The  cross  section  for  ^B(p,n)^C  reaction  has  been  studied  in 


the  past. 


Discrepancies  by  as  much  as  a  factor  of  four  exist  in 


published  total  cross  sections.  Thin  target  techniques,  involving  the 
evaporation  of  a  thin  layer  of  boron  on  a  surface  have  been  extensively 
used  in  the  past.  The  discrepancies  in  past  measurements  were  probably 
caused  by  uncertainties  in  composition  or  thickness  of  these  thin 
targets.  We  have  tried  to  avoid  the  first  problem  by  using  thick  tar- 


The  ^B(p,n)^C  reaction  is  an  interesting  candidate  for  boron 


depth  profiling  experiments.  The  cross  section  for  production  of  neutrons 


3 


is  relatively  large,  in  the  ten  millibarn/steradian  region,  and  varies 
slowly  with  proton  energy.  Only  ground  state  neutrons  can  be  produced 
at  proton  energies  below  5.2  MeV.  The  reaction  leads  to  reasonable 
sensitivity  and  resolution  for  profiling,  as  discussed  in  detail  in 
Chapter  VI.  From  an  applied  physics  point  of  view,  boron  profiling 
could  be  used  by  the  semiconductor  industry,  where  boron  is  used  in 
making  p-type  materials  by  implanting  specific  concentrations  of  boron 
at  specific  depths  in  intrinsic  or  n-type  materials. 

The  first  objective  was  accomplished  by  studying  neutron  yields 
from  thick  and  B  metal  targets  through  neutron  Time- of- Flight 

(TOF)  techniques,  after  carefully  determining  the  neutron  detector 
efficiency.  After  resolving  the  cross  section  discrepancies,  the  second 
objective  was  met  by  profiling  known  amounts  of  boron  implanted  in 
silicon  at  known  depths.  The  feasibility  of  the  reaction  technique  for 
profiling  was  then  determined  by  comparing  the  experimental  resolution 
and  sensitivity  to  theoretical  expectations. 


4 


II.  EXPERIMENTAL  METHODS 

A.  Thick  Target  Kinematics. 

Thick  boron  targets  of  about  1  mm  thickness  were  bombarded  with 
a  nanosecond-pulsed  proton  beam.  The  thickness  of  the  target  is 
important,  because  the  proton  beam,  after  striking  the  target  with 
an  incident  energy,  must  be  degraded  in  energy  to  below  the  reaction 
threshold  prior  to  exiting  the  target. 

The  main  advantage  of  the  thick  target  technique  is  that  the 
cross  section  for  the  production  of  neutrons  as  a  function  of  proton 
energy  (E^)  can  be  obtained  from  only  the  one  incident  proton  beam 
energy.  This  technique  also  has  the  advantage  that  targets  made  from 
bulk  material  should  be  more  easily  reproduced  than  chin  targets. 

A  proton  beam  incident  on  a  thick  boron  target  will  interact  and 
slow  down  in  accordance  with  the  stopping  power  of  the  target  material. 
Therefore,  a  continuum  of  proton  energies  from  incident  to  below  thres¬ 
hold  is  produced.  Because  there  is  a  possibility  of  a  nuclear  reaction 
at  all  energies  above  threshold  there  will  also  be  a  continuum  of  neutron 
energies  produced.  Considering  the  ^B(p,n)^'*'C  nuclear  reaction  as  a 
two  body  process,  conservation  of  energy  and  of  momentum  give  rise  to 
"kinematic  restrictions"  of  the  reaction.  The  resulting  energy  of  the 
neutron  (E^)  is  a  function  of  its  emission  angle  (0),  energy  of  Che 
proton  and  masses  of  the  four  interacting  particles. 

As  a  result  of  the  negative  Q  value  the  ^B(p,n)^C  nuclear  reaction 
has  the  following  kinematic  restrictions.  The  reaction  has  a  kinematic 
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threshold  value  of  E  *  3.016  MeV  where  neutrons  are  emitted  with 

P 

E  “20  KeV  only  at  0,  .  *  0°.  As  E  increases  neutrons  are  emitted  at 
n  lab  p 

greater  angles  until  finally  a  "cone  opening"  energy  of  3.042  MeV  is 

reached.  Neutrons  are  then  emitted  at  all  angles  in  the  lab  frame. 

Values  for  E  from  threshold  to  cone  opening  can  give  rise  to  double 
P 

values  of  E  in  the  lab  frame  corresponding  to  forward-  and  backward- 
n 

emitted  neutrons  in  the  center  of  mass  frame.  The  backward-emitted 

neutrons  will  have  a  lab  E^  *  20  KeV  at  threshold  and  zero  energy  at 

the  cone  opening  energy.  At  Ep  greater  than  cone  opening  energy  only 

one  E  is  allowed  for  a  given  0,  ,  and  E  .  For  example,  at  E  * 
n  6  lab  p  p 

4.5  MeV  and  0.  »  0°,  E  =  1.69  MeV  while  at  the  same  E  with 

lab  n  p 

9.  ,  =  180°,  E  =  0.87  MeV. 
lab  n 


B.  Transformation  Equations. 


A  neutron  TOF  spectrometer  was  used  and  the  number  of  neutrons 
detected  as  a  function  of  their  flight  time  t  from  the  target  was 
measured.  More  precisely,  the  number  of  neutrons  detected  N^t(t)At  is 
the  number  of  neutrons  detected  at  time  t,  summed  over  time  interval 
At.  Knowing  the  experimental  time  of  flight  and  the  flight  path  from 
the  target  to  the  detector  the  corresponding  neutron  energy  can  be  cal¬ 
culated.  Nonrelativistically  the  kinetic  energy  of  the  neutron  is 
2  2 

^  mX  /  t  ,  where  X  is  the  flight  path  and  m  is  the  mass  of  the  neutron. 
The  neutron  energy  interval  AE^  which  corresponds  to  At  is  related  to  a 
proton  energy  interval  AE^  through  the  kinematic  equations. 
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t 


I 


» 


k  continuous  neutron  TOF  spectrum  can  be  transformed  to  a  neutron 
11 

energy  spectrum. 


N(En,9)  - 


Ndet(C)At 


2.1 


N  •  AE 
P  n 


AG  •  e(E  ) 
n 


where  N(En>9)  is  the  number  of  neutrons  emitted  at  angle  0  per  incident 
proton,  unit  neutron  energy  interval  and  solid  angle.  The  neutron 


detector  subtends  solid  angle  AG  at  the  target,  is  the  number  of 

incident  protons  that  strike  the  target  and  e(E  )  is  the  neutron  detector 

n 


efficiency. 

The  detector  efficiency  is  defined  as  the  number  of  neutrons 

detected  divided  by  the  number  of  neutrons  that  strike  the  detector.  In 

present  work,  e(E  )  is  experimentally  determined  using  the  ^Li(p,n)^Be 
n 

nuclear  reaction  as  the  standard.  The  experimental  neutron  detector 
efficiency  is  discussed  in  detail  in  Chapter  IV. 

Once  eCE^)  Is  experimentally  determined  the  differential  cross 
section  a(0)  can  be  calculated  from11 


o  (0) 


,1  dEj 

^n  dx 


Ndet(t)« 


N  •  AE  •  AG  •  e(E  ) 
p  p  n 


2.2 


i 
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•  (-  j=9-)  •  M  +  a(-  • 

pdx  d _  d  p  dx  a 


(0.S022) 


where  "b"  denotes  the  number  of  atoms  of  boron  (B)  in  the  target 

molecule  and  "a"  is  the  number  of  atoms  of  the  other  element  (A)  in  the 

molecule.  For  example,  a  target  would  have  b  -  2  and  a  *  3.  The 

1  dE 

stopping  power  of  the  given  element  is  given  by  —  -j^P  as  compiled  by 
Janni, M  is  the  corresponding  element's  molecular  weight  and  T  is 
Avogadro's  number.  The  number  0.8022  accounts  for  the  fact  that  boron 
is  80.22%  UB  and  19.88%  10B. 

The  stopping  cross  section  per  22B  atom  must  be  calculated  at 

some  energies  because  Janni  only  gives  the  stopping  power  for  each 

element  at  200  KeV  intervals.  The  interpolation  formula  used  is  a 

14 

generalized  form  of  the  Bethe-Bloch  equation, 

—  4^P  *  C  •  [ln(E  )  +  C,]  /E  .  2.4 

ndx  1  p  2  p 

The  constants  and  C,  were  obtained  by  solving  the  above  equation 

1  dE 

simultaneously  using  values  of  —  jjjP  from  Eq.  2.3  at  proton 

1  dE 

energies  of  3.0  and  4.6  MeV.  Values  of  — t=P  at  intermediate  energies 

n  dx 

agree  within  0.1%  with  the  values  obtained  from  Eq.  2.3  and  Janni' s 

tabulated  stopping  powers  over  the  proton  energy  range  from  3.0  to  4.6 

-21 

MeV.  For  B?0^,  values  of  constants  are  -  5.525  x  10 
.  2  2 , 

MeV  cm  /ll^  atom*  ^2  *  while  for  elemental  boron, 

C.  -  1.644  x  10-21  MeV2cm2/ll  and  C,  -  2.982. 

i.  B  atom  2 


III.  EXPERIMENTAL  DETAIL 


A.  Hardware. 

The  same  experimental  set  up  as  described  by  Burke^  was  used  in 
the  determination  of  the  ^B(p,n)^C  cross  section.  A  more  detailed 
discussion  of  the  following  areas  can  be  found  in  his  thesis:  pulsed 
beam  mechanism,  target  holder,  target  chamber,  neutron  shielding, 
neutron  detector,  and  electronic  set  up. 

The  University  of  Oregon  Van  de  Graaff  accelerator  was  used  to 
produce  a  4.6  MeV  pulsed,  klystron  bunched,  proton  beam.  Deflection 
and  focussing  magnets  were  used  to  direct  the  beam  through  three 
collimators  to  the  target  on  the  32°  East  beam  line.  The  resulting 
beam  size  was  approximately  3  mm  in  diameter  and  had  a  beam  pulse 
repetition  period  of  4  to  16  psec  depending  on  the  type  of  target  used. 
An  overall  time  resolution  of  1.5  nsec  was  obtained.  Time  average  beam 
currents  were  about  0.5  microamperes.  The  beam  line  vacuum  near  the 
target  was  ^5  x  10  ^  torr. 

Two  types  of  targets  were  used  in  the  determination  of  the  cross 

i k 

section:  ®2^3  anc*  B  netal.  Analytical  reagent  grade  Boric  Acid  was 
used  to  make  the  targets,  through  the  following  reaction, 

2  H,B0,  7— ►  Bo0,  +  3  H,0  . 

3  3  4  2  3  2 

* 

Mallinckrodt  Chemical  works,  Boric  Acid  granular,  Heavy  metal  to 
0.01%  max  impurity. 


9 


Boric  acid  was  heated  in  a  stainless  steel  target  holder  to  a  red 
heat  in  a  partial  vacuum.  As  the  sample  was  heated,  the  pressure 
inside  the  bell  jar  was  decreased  slowly  to  a.25  microns.  Extreme  care 
had  to  be  taken  to  insure  most  of  the  water  was  driven  off,  resulting 

in  a  clear,  colorless,  glass-like  solid.  Boron  metal  was  commercially 

* 

purchased  from  Alfa  Ventron  as  a  polycrvstalline  solid.  Because 
boron  metal  has  a  hardness  rating  of  9.5  (MOHS)1^  a  diamond  toothed 
saw  had  to  be  used  to  cut  the  boron  chunks  into  target  thicknesses 
of  1  mm.  Two  different  types  of  targets  were  prepared.  One  target  had 
the  front  surface  untouched,  the  other  was  sanded  smooth  using  a  diamond 
impregnated  sanding  disk.  Both  targets  were  mounted  on  tantalum  target 
holders  using  epoxy  glue. 

Two  different  types  of  target  holders  were  used.  The  closed  end 
stainless  steel  "top  hat"  used  in  the  preparation  of  the  B,0^  targets, 
had  an  inside  diameter  of  1  cm,  wall  thickness  of  0.25  mm  and  an  overall 
length  of  5  cm.  The  boron  metal  targets  were  mounted  in  tantalum 
dishes,  0.6a  mm  thick  and  1.4  cm  inside  diameter.  Both  target  holders 
were  hermetically  sealed  to  the  evacuated  target  chamber  by  an  "o-ring" 
and  held  in  place  by  atmospheric  pressure. 

The  target  chamber  was  a  stainless  steel  tube  with  an  inside 
diameter  of  1.27  cm  and  an  overall  length  of  lb. 5  cm.  The  chamber 
was  electrically  isolated  from  the  target  holder  by  a  1  cm  long  glass 
tube  and  also  from  the  grounded  beam  pipe  by  a  non-conducting  pipe- 


Alfa  Division,  Ventron  Corp.  Boron  pieces  crystalline,  5  cmx  down 
99.3"  pure  with  respect  to  metallic  impurities. 
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coupling  flange.  The  target  chamber  was  biased  at  -300  volts  with 
respect  to  the  target  holder  for  secondary  electron  suppression. 

The  neutron  shielding  which  surrounded  the  target  was  made  by 
Lunnon^  and  designed  to  reduce  room  background.  The  system  consisted 
of  two  metal  tanks,  1  meter  in  diameter,  containing  a  lithium  carbonate  - 
water  slurry  and  nine  paraffin  -  lithium  carbonate  wedge-shaped  sections. 
The  pie  sections  were  mounted  between  the  two  tanks  on  a  horizontal 
plane  parallel  to  the  target  chamber.  One  of  the  pie  sections  was 
made  of  high  density  polyethylene  and  had  a  horizontal  conical  hole 
cut  in  its  center  with  half  angle  of  1/2°,  allowing  neutrons  emitted 
from  the  target  to  strike  the  detector  unobstructed. 

The  neutron  detector  consisted  of  a  cylindrical  Naton  136  plastic 
scintillator,  1.27  cm  thick  and  6.2  cm  in  diameter,  optically  coupled 
to  a  RCA  C70133  B  photomultiplier.  The  neutrons  were  detected  through 
the  hydrogen  and  carbon  atom  recoils  in  the  plastic  scintillator.  The 
entire  assembly  was  magnetically  shielded  to  decrease  the  angular 
dependence  of  the  photo  tube  gain  to  less  than  the  1 %  level. 

A  diagram  of  the  electronics  used  in  the  measurement  is  shown  in 
Fig.  2.  The  analog  to  digital  conversion  of  the  Time-to-Ampiitude 
Converter  (TAC)  output  pulses  was  done  by  a  dual  twelve  bit  Nuclear 
Data  161F  Analog  to  Digital  Converter  (ADC). 
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B.  Data  Acquisition. 

The  pulse-height  electronics  were  standardized  to  a  pulse  height 
241 

spectrum  of  Am  >  rays  prior  to  each  day's  run.  The  electronics  were 

checked  several  times  during  the  day's  run  and  corrected  for  any  drift. 

241 

The  detection  (or  "slow")  bias  was  set  to  1/17  of  the  “  Am  60  KeV 
Y  ray  full  energy  peak.  The  timing  (or  "fast")  bias  was  set  below 
this  level.  The  fast  bias  was  always  adjusted  such  that  a  pulse  height 
spectrum  gated  with  only  this  bias  matched  the  spectrum  gated  only  by 
the  slow  bias  for  pulse  heights  above  the  level  of  the  slow  bias. 

Neutron  TOF  spectra  were  taken  in  three  groups  over  a  one  year 
period.  The  B,0^  targets  were  used  to  measure  the  neutron  yields  at 
the  following  eleven  angles:  0°,  20°,  31°,  50°,  69.8°,  80.3°,  90°, 

100°,  114°,  138°  and  155°.  The  flight  path  varied  depending  on  the  room 
geometry.  The  TOF  spectra  at  0°  to  31°  were  taken  at  2.19  meters  while 
those  taken  at  50°  to  155°  were  at  1.91  meters.  Boron  metal  targets 
were  also  used  to  measure  TOF  spectra  at  0°.  No  appreciable  difference 
was  noticed  between  the  two  types  of  B  metal  targets.  The  ADC  TOF 
processing  dead  time  was  v 8Z  for  B^0^  targets  and  \-lb2  for  B  metal 
targets.  This  dead  time  effect  was  compensated  for  by  using  the  ex¬ 
ternal  clock  input  port  of  the  ADC,  which  timed  each  run  for  a  pre¬ 
selected  amount  of  beam  charge.  The  number  of  incident  protons  that 
struck  the  target  was  derived  from  the  resulting  total  integrated  charge. 
This  charge  was  an  experimental  variable  and  was  chosen  to  give  counting 
statistics  of  •vlX  with  ^10&  counts  above  the  fast  bias.  The  integrated 
charge  for  the  B.,0^  targets  was  200  uc  at  each  angle  and  for  the  B 
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metal  targets  It  was  100  uc. 

The  resulting  TOF  spectra  were  stored  in  a  1024  channel  array 
of  number  of  counts  versus  flight  time  in  a  PDP-7  computer.  After 
completion  of  each  run  the  spectrum  was  written  on  magnetic  DEC  tape. 

A  TOF  spectrum  for  the  *^B(p,n)^C  reaction  with  a  B.,0^  target, 
with  an  incident  proton  energy  equal  to  4.517  MeV  and  laboratory 
angle  equal  to  zero  degrees,  is  shown  in  Fig.  3.  The  flight  time  is 
measured  from  right  to  left  with  channel  880.5  corresponding  to  zero 
flight  time.  The  prompt  y  ray  peak  is  located  in  channel  873.  The 
sharp  leading  edge  of  the  neutron  production  is  located  in  channel  755 
and  corresponds  to  the  maximum  neutron  energy  of  1.69  MeV.  The  neutron 
energy  corresponding  to  a  given  channel  number  C  is 

E  _  _ I2  3  1 

n  L880.5-CJ  ’  3,1 

where  the  constant  163  incorporates  a  flight  path  of  2.197  meters  and 
a  channel  calibration  of  0.9756  nsec  per  channel. 

In  contrast  to  previous  TOF  spectrum,  a  TOF  spectrum  from  the  same 

target  taken  with  an  incident  proton  beam  below  the  threshold  energy 

revealed  a  completely  flat  background  with  a  similar  prompt  y  ray  peak, 

and  the  same  smaller  peak  at  channel  810.  The  smaller  peak  is  thought 

to  be  a  result  of  a  higher  beam  mass  or  an  unlikely  gating  difficulty 

visible  only  when  y  rays  are  intense.  The  flat  background  indicates 

that  no  reactions  with  higher  Q  values  contaminate  the  spectrum.  Back- 

IS  L3 

ground  effects  from  the  0(p,n)  F  reaction  were  Investigated  at 


— 


■ 
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Ep  ■  4.517  MeV,  0^^  “  0°  withanSiO,  target.  The  resulting  TOF  spec¬ 
trum  is  shown  in  Fig.  4.  Again  the  same  two  y  ray  peaks  are  seen  in 
channels  beyond  800.  Some  structure  is  seen  in  the  region  from 
channels  650  to  800,  with  a  maximum  value  of  70  counts  and  a  back¬ 
ground  level  of  v30  counts.  Comparing  this  spectrum  to  the  B.,0^  TOF 
spectrum  one  can  see  that  the  statistical  uncertainty  in  the  number  of 

counts  from  the  B^0^  target  is  approximately  equal  to  the  number  of 

IS 

counts  from  the  SiO.,  target.  Thus  the  effect  of  0  in  the  B.,0_ 

i.  -  j 

target  is  very  small  and  is  neglected  in  this  experiment. 

A  total  of  fifteen  \i(p,n)^Be  neutron  TOF  spectra  were  taken  on 
two  days,  three  months  apart.  These  TOF  spectra  were  used  to  determine 
the  neutron  detector  efficiency,  which  is  discussed  in  detail  in  the 
next  chapter.  Lithium  metal  targets  were  prepared  as  described  by 
Burke'-*'  and  mounted  in  the  top  hat  target  holder.  The  same  electronic 
set  up  was  used  with  a  flight  path  of  2.14  meters,  total  integrated 
charge  of  25  pc  and  an  ADC  dead  time  of  20*. 

The  time  per  channel  calibration  was  taken  each  day  upon  completion 
of  the  day's  run,  using  an  Ortec  462  Time  Calibrator.  Several  spectra 
were  taken  at  various  dead  times  corresponding  to  the  experimental 
dead  times  noted  during  the  B^0^,  B  metal  or  Li  targets  run  on  that 
day.  The  average  channel  width  was  M3.97  nsec  per  channel. 
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IV.  DATA  ANALYSIS  FOR  CROSS  SECTION 


A.  Determination  of  Neutron  Detector  Efficiency. 


The  neutron  detector  efficiency,  as  defined  in  Chapter  II,  can 
be  expressed  by  rewriting  Eq.  2.1, 


“E„> 


N(E  ,0)  •  AE  •  Aft  •  N 
n  n  p 


The  neutron  detector  efficiency  has  been  studied  by  Burke 
utilizing  the  same  electronic  set  up.  The  detector  efficiency  is  sub¬ 
ject  to  experimental  conditions  and  as  such  should  be  experimentally 
determined.  For  example,  the  efficiency  of  the  detector  is  a  function 
of  the  signal  processing  bias  settings  (fast  and  slow),  the  quality  of 
the  optical  coupling  between  the  scintillator  and  the  photomultiplier 
and  sometimes  the  age  of  the  scintillator. 

The  \i(p,n)'Be  nuclear  reaction  is  used  to  determine  e(E  )  bv 

n 

experimentally  determining  N^^,  AE^,  Aft  and  N^.  The  number  of 

neutrons  emitted  in  the  ^Li(p,n)^Be  reaction,  N(E  ,6),  can  be  calculated 

n 

if  the  (p,n)  cross  sections  and  atomic  stopping  cross  sections  are 


known.  The  number  of  neutrons  emitted  can  be  solved  from  Eq.  2.1 
and  2.2, 


N(E  ,0)  -  q _  dEp  4.2 

(1/n) (dE  /dx)  dE 
P  n 

where  the  transformation  derivative  dE  /dE  is  determined  through  the 

P  n 

kinematic  relationships  for  the  \i(p,n)^Be  reaction  and  j-j-  is  the 
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atomic  stopping  cross  section  per  7Li  atom.  The  case  of  the  7Li(p,n)7Be 

reaction  is  more  complicated  chan  implied  by  Eq.  4.2  because  two 

groups  of  neutrons  can  be  emitted  if  E  is  above  2.378  MeV. 

P 

The  cross  section  for  the  Li(p,n)7Be  reaction  at  0.  ,  *  0° 


lab 


11,17-19 


has  been  studied  extensively  and  has  been  reported  in  detail. 

Lefevre^  has  used  Eq.  4.2,  the  ground  state  cross  sections  of  Burke^ 


and  the  first  excited  state  cross  sections  inferred  from  the  measure- 
17  18 

ments  of  Bevington  and  Borchers  to  calculate  N^CE^.00)  for  ground 

state  neutrons  and  N.  (E  ,0°)  for  first  excited  state  neutrons. 

1  n 

The  number  of  detected  neutrons  in  a  time-of-flight  experi¬ 


ment  is  converted  to  a  modified  energy  spectrum 

N  (E  ,0°)  -  -  . 

exp  '  n’  '  N  •  AE  •  A il  ' 
p  n 


4.3 


This  is  just  Eq.  2.1  with  the  efficiency  efE^) 
efficiency  can  now  be  determined  from 


1.  The  detector 


c<En) 


N  (E  ,0") 
exp  n 


N  (E  ,0°)  +  N  (E  ,0°) 
on  in 


4 . 4 


Lithium  TOF  spectra  for  incident  proton  energies  from  2.9  MeV  to  5.1 

MeV  were  transformed  to  Ngxp(En,0°)  spectra  using  the  program  T0F-NE, 

written  by  Burke^  (Appendix  C).  Utilizing  Eq.  4.4,  an  efficiency 

table  was  calculated  for  each  transformed  Li  TOF  spectra.  The  number 

of  neutrons  emitted  N  (E  ,0°)  +  N, (E  ,0°)  was  calculated  separately 

on  in 

at  each  incident  energy  from  Lefevre's  tables.  Efficiencies  were 
calculated  for  each  of  the  two  separate  days  and  averaged,  taking  into 
account  the  target  front  surface  contamination  effects.  Average 
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efficiencies  for  each  day  were  within  the  experimental  statistical 
errors,  but  the  two  days  differed  by  ±5%  from  the  overall  average 
or  by  somewhat  more  than  the  statistical  uncertainties. 

Extreme  care  must  be  exercised  in  the  steering  of  the  pulsed 
proton  beam  onto  the  target.  Once  the  beam  is  properly  steered 
onto  the  center  of  the  target,  over  or  under  focusing  of  the  beam 
will  have  little  effect  on  the  resulting  TOF  spectrum.  The  prepara¬ 
tion  of  the  Li  is  also  very  important.  Target  front  surface 
contamination  effects  extending  as  far  as  300  KeV  in  proton  energy 
loss  were  noted  in  efficiencies  determined  for  3.0  and  5.1  MeV  protons 
incident  on  the  same  target. 

The  resulting  overall  efficiency  for  the  neutron  detector  and 

that  calculated  by  Burke^  are  shown  in  Fig.  5.  The  two  sets  of 

efficiencies  agree  to  ±4%  for  greater  than  0.4  MeV.  For  E^  less 

than  this,  Burke's  values  are  lower  by  more  than  7%.  In  the  neutron 

energy  interval  of  0.2  MeV  to  0.4  MeV  the  ^Li(p,n)^Be  cross  section 

at  0,  .  ■  0°  has  a  very  sharp  minimum.  It  is  possible,  in  this  region, 
lab 

that  neutron  scattering  could  lead  to  significant  effects  in  the 
efficiency  calculation.  The  dip  in  the  efficiency  table  at  250  Kev 
is  due  to  Li  resonance  scattering  at  that  energy.  Burke's  results 
were  corrected  for  this  effect. 

The  efficiency  table  used  in  this  work  was  based  on  Burke' s^ 
measured  cross  section  values.  Because  Burke  used  the  same  experimental 
set  up  it  is  hoped  that  experimental  systematic  errors  will  tend  to 
cancel.  A  comparison  of  Burke's  published  zero  degree  ground  state 


50 
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cross  sections  and  those  tabulated  by  Liskien  and  Paulsen  are  shown 
in  Fig.  6.  Burke's  values  are  not  included  in  those  reported  by 
Liskien  and  Paulsen  and  are  ^10%  higher. 

The  neutron  detector  efficiency  used  through  the  remainder  of 
this  work  has  an  estimated  overall  accuracy  of  i7Z.  If  the  cross 
section  values  by  Liskien  and  Paulsen  were  used  the  efficiencies 
would  be  larger  by  M.0%.  Thus,  there  is  a  possible  systematic  error 
of  ■v-lOZ  in  the  efficiencies  used. 

B.  Conversion  of  Time-of-Flight  Spectra  to  Cross  Section. 

A  constant  flat  background  was  subtracted  from  each  ^B(p,n) 
neutron  TOF  spectrum.  Each  spectrum  was  then  smoothed  to  increase 
counting  statistics  by  using  a  three  channel  running  average.  The  TOF 
spectrum  was  then  labeled  in  channel  0-5  with  experimental  parameters 
that  are  used  in  the  transformation  computer  programs. 

The  prepared  TOF  spectra  were  transformed  to  laboratory  differen¬ 
tial  cross  sections  versus  proton  energy  by  the  program  B11PNX  (Appendix 
D) .  The  computer  program  was  a  modified  version  of  LITPN1  written  by 
Burke.  The  program  uses  Eq.  2.2  and  efficiencies  described  in  the 
previous  section.  Time-of-f light  spectra  for  both  B,0,  and  B  metal 

i  i 

targets  were  transformed  using  the  B11PNX  program  where,  of  course, 
appropriate  atomic  stopping  cross  sections  were  used  for  each  target. 


a 


neutron  production  cross  sections. 
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V.  CROSS  SECTION  ANALYSIS 

A.  Results  and  Discussion. 

The  ^B(p,n)^C  zero  degree  differential  cross  section  for  B.,0, 

and  3  metal  targets  are  compared  to  the  published  results  by  Van  der 
9 

Zwan  and  Gieger  in  Fig.  7.  The  target  was  bombarded  with  .  5 1 7 

MeV  incident  protons,  while  4.618  MeV  protons  were  used  for  the  3 
metal  target.  The  disagreement  between  the  two  boron  targets  is  less  than 
±4*  but  is  larger  than  uncertainties  due  to  counting  statistics.  The 
fact  that  the  B-,0^  target's  cross  section  was  lower  at  high  energies 
and  larger  near  4.1  and  3.5  MeV  could  be  due  to  two  types  of  target 
impurities.  The  B.,0^  target  could  have  absorbed  water  molecules  on  its 
front  surface  during  transportation  of  the  target  from  the  bell  jar  to 
the  target  chamber.  The  B  metal  target,  a  polycrvstalline  solid,  could 
have  had  slight  impurities  on  the  surface  on  internal  crystals.  There¬ 
fore,  as  the  proton  beam  penetrated  further  into  the  target  it  could 
have  passed  through  pockets  of  impurities,  thereby  giving  lower  cross 

section  results.  Both  the  3^0 ^  and  B  metal  target  gave  cross  section 

9 

values  larger  than  those  reported  by  Van  der  Zwan  and  Gieger.  .An 
experimental  error  bar  for  the  latter's  work  is  shown  in  Fig.  7  at 
3.58  MeV. 

The  ^3(p,n)^C  differential  cross  sections  measured  at  lab  angles 
of  0°,  20°,  31°,  50°,  69.8°,  30.3°,  90°,  100°,  114°,  138.2°,  and  155° 
are  shown  in  Fig.  8.  These  cross  sections  were  obtained  from  bom¬ 
barding  a  B,0^  target  with  a  4.517  MeV  incident  proton  beam.  The 


E 


P,  Lab 


(  M  e  V  ) 


B(p,n)  C  differential  cross  sections  for  laboratory 
angles  of  0  through  155°.  The  reaction  angle  of  each  curve 
is  labeled  on  the  right.  Each  curve  is  labeled  on  the  left 
side  with  a  vertical  offset  (in  parentheses)  to  be  subtracted 
from  the  ordinate  to  obtain  the  cross  section  at  that  angle. 
The  reaction  threshold  is  at  3.015  MeV. 
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numerical  values  of  the  cross  section  and  their  statistical  uncertain¬ 
ties  are  given  in  Appendix  A.  The  pronounced  dip  in  the  cross  section 
at  3.33  MeV  at  0°,  increasing  in  proton  energy  to  3.62  MeV  at  100°, 
corresponds  to  a  neutron  energy  of  0.433  MeV.  Because  this  dip  was 

always  at  the  same  neutron  energy,  it  is  probably  due  to  resonance 

16  2 1 

scattering  by  the  target  0  nuclei.  A  closer  look  at  results 

from  the  B^O^  and  B  metal  targets  in  Fig.  7  shows  that  the  B  metal 
target  values  do  not  exhibit  such  a  drastic  dip.  A  comparison  of  the 
eleven  curves  in  Fig.  8  shows  that  the  zero  degree  cross  section  is 
the  largest. 

Angular  distributions  in  the  center  of  mass  system  were  fitted 
with  a  series  of  Legendre  polynomials; 

N 

a  (9)  =  Z  \  P  (cose  )  . 

L=0 

The  conversion  to  the  center  of  mass  system  and  calculation  of  the 

Legendre  polynomial  coefficients  was  done  using  the  program  LSQFT3 

(Appendix  E) .  This  program  was  a  modified  version  of  Burke's^  LTHM3. 

The  LSQFT3  program  allows  the  operator  to  input  at  the  teletype  computer 

interface  the  masses  and  Q  value  of  the  two  body  reaction  under  study, 

the  number  of  Legendre  polynomials  coefficients  to  be  calculated  (up 

to  6)  and  the  number  of  differential  laboratory  cross  sections  to  be 

used  as  a  basis  for  calculations  (up  to  15). 

The  total  cross  section,  4TrAQ,  for  the  ^B(p,n) reaction  is 

9 

compared  to  results  published  by  Van  der  Zwan  and  Geiger  and  to 

4 

normalized  results  of  Gibbons  and  Macklin  in  Fig.  9.  Because  the 


«u*d  »#„  —  S  Z  I 


present  results  with  error  bars  representing  the  statistical  uncertainties 
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latter  authors  were  unable  to  determine  their  target  thickness,  their 
published  results  were  normalized  to  results  published  by  Blaser  et  al.“ 
Thus  the  results  of  Gibbons  and  Macklin  were  multiplied  by  a  factor 
of  1.99  to  normalize  them  to  the  present  work  near  the  3.6  MeV  peak. 
Statistical  uncertainties  for  the  present  work  are  shown  with  error 
bars.  Van  der  Zwan  and  Geiger's  results  are  plotted  at  100  KeV  intervals 
and  show  a  peak  value  at  3.58  MeV.  Their  experimental  uncertainty  is 
shown  by  an  error  bar  at  this  energy.  Within  the  uncertainties,  there 
is  agreement  between  Van  der  Zwan  and  Geiger's  results  and  the  present 
work.  The  numerical  values  of  the  total  cross  section  and  the  associated 
uncertainties  are  given  in  Appendix  B. 

The  Legendre  polynomial  power  series  coefficients  A^  through  A^ 

for  the  present  work  and  the  values  reported  by  Van  der  Zwan  and 
9 

Geiger  are  shown  in  Fig.  10.  The  coefficients  and  associated  statis¬ 
tical  errors  for  the  present  work  are  tabulated  in  Appendix  B.  The 
error  bars  in  Fig.  10  represent  the  statistical  uncertainty  for  the 
present  work.  Van  der  Zwan  and  Geiger’s  values  are  plotted  at  100  KeV 
intervals  starting  at  3.5  MeV.  They  report  an  isotropic  behavior  from 
threshold  to  3.5  MeV,  but  five  polynomials  were  required  to  fit  their 
data  between  3.5  and  4.5  MeV.  The  present  results  agree  with  a  fit  of 
five  polynomials  at  higher  energies  but  do  not  show  the  isotropic 
behavior  at  low  energies.  Increasing  the  number  of  polynomials  to  six 
did  not  improve  the  quality  of  the  fit.  Due  to  the  kinematic  restric¬ 
tions  in  the  ^B(p,n)^C  reaction,  the  coefficients  calculated  below 
3.2  MeV  are  not  reliable  because  the  back  angle  data  does  not  exist 


(mb/sr)  Ap  (mb/sr)  Ai  (mb/sr)  A4  (mb/sr) 
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or  there  are  large  uncertainties  due  to  the  low  neutron  detector 
efficiency.  At  values  greater  than  3.8  MeV  the  present  data  show  a 
general  agreement  with  those  of  Van  der  Zwan  and  Geiger. 

The  center-of-mass  cross  sections  as  a  function  of  center-of-mass 
angle  for  selected  laboratory  proton  energies  are  shown  in  Fig.  11. 

The  solid  curves  represent  calculated  cross  sections  based  on  the  five 
polynomial  coefficients  shown  in  Figs.  9  and  10.  The  cross  sections 
appears  to  be  forward  peaked  at  higher  proton  energies. 

B.  Uncertainties. 

The  uncertainties  in  this  experiment  can  be  divided  into  two  areas: 
those  associated  with  the  cross  section  and  those  associated  with  the 
proton  energy. 

The  uncertainty  in  the  cross  section  can  be  determined  from 
Eq.  2.2  by  taking  the  square  root  of  the  sum  of  the  square  of  the  un¬ 
certainties  in  each  of  the  terms.  The  number  of  neutrons  detected 
was  subject  to  statistical  uncertainties  in  the  number  of  counts  (^1% — 
which  is  a  function  of  the  flight  time)  and  in  the  time  correlated 
background  (M.%) .  The  resulting  uncertainty  for  N(jet  was  about  1.4%. 

The  remainder  of  the  uncertainties  in  the  other  terms  were  due  to 
systematic  errors.  The  stopping  cross  section  had  an  overall  uncertainty 

of  6%.  This  value  was  due  in  part  to  Janni's^  elemental  — 

r  0  dx 

("v4 %)  and  possible  target  impurities  that  could  give  rise  to  another 
4%  error.  The  efficiency  of  the  neutron  detector  has  an  associated 
uncertainty  of  7%,  although  the  uncertainty  is  probably  lower  at  higher 
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Fig.  11.  ^B(p,n)^C  angular  distributions  of  the  center  of  mass  system. 

The  laboratory  proton  energy  of  each  curve  is  labeled  on  the 
left.  The  top  four  curves  here  are  labeled  on  the  right  side 
with  a  vertical  offset  (in  parentheses)  to  be  subtracted  from 
the  ordinate  to  obtain  the  cross  section  at  the  given  energy. 
The  solid  curves  represent  the  calculated  cross  section  as 
determined  from  the  Legendre  polynomials. 
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energies  as  discussed  in  Chapter  IV.  The  number  of  incident  protons 
had  an  associated  uncertainty  of  0.3%  due  in  part  to  beam  current  in¬ 
tegration  uncertainty.  The  solid  angle  uncertainty  was  due  to  the 
combined  effects  of  the  area  of  the  scintillator  and  the  flight  path 

which  gave  a  0.9%  uncertaintv.  Finally  the  uncertainty  in  dE  was 

•  n 

^2%  due  mainly  to  the  channel  width  uncertainty. 

Combining  the  uncertainties  as  mentioned  above,  a  9.6%  uncertainty 

would  be  associated  with  the  determination  of  the  cross  section.  If 

the  neutron  efficiency  were  measured  with  respect  to  Liskien  and 
19 

Paulsen's  data,  as  discussed  in  Chapter  IV,  results  would  be  con¬ 
sistently  lower  by  about  10%. 

Uncertainties  associated  with  the  proton  energy  axis  of  the  cross 

section  curves  depend  on  flight  path  and  flight  time  uncertainties, 

and  kinematic  transformation  of  the  flight  time  to  proton  energy. 

The  flight  path  had  an  uncertainty  of  0.6%  due  mainly  to  the  thickness 

of  the  scintillator.  The  flight  time  uncertainty  was  a  function  of  the 

overall  timing  resolution  (FWHM  *  1.5  nsec)  and  the  Ortec  time  mark 

generator  which  gave  the  number  of  nsec  per  channel.  At  neutron  energies 

greater  than  the  bias  settings,  -  1,  which  meant  that  the  uncertainty 

P 

for  both  the  neutron  and  proton  energy  were  about  the  same.  Taking 
all  this  into  account,  uncertainties  in  proton  energy  at  3.2  and  4.5 
MeV  were  ±3  KeV  and  ±21  KeV  respectively. 
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C.  Summary. 


The  first  objective  of  this  thesis  has  been  satisfactorily  met. 

Table  I  is  a  modified  summary  of  past  work  as  first  described  by  Van 

9 

der  Zwan  and  Geiger.  There  is  excellent  agreement,  within  experi¬ 
mental  errors,  between  Van  der  Zwan  and  Geiger  and  the  present  work 
for  the  total  cross  section.  This  was  very  significant,  because  both 
the  thin  target  (Van  der  Zwan  and  Geiger)  and  thick  targets  (present 
work)  techniques  gave  similar  results.  Overall  the  present  work 
agrees  well  with  all  their  results  except  at  zero  degrees,  where  there 
was  a  26%  difference.  The  reason  for  this  discrepancy  is  not  known. 
However,  Marion  et  al.^  zero  degree  cross  section  at  *  3.46  MeV, 
which  was  obtained  using  thin  target  techniques,  was  in  excellent 
agreement  with  present  results. 

Careful  study  of  zero  degree  cross  sections  from  and  B  metal 

targets  revealed  that  there  is  excellent  agreement  between  the  two  at 
low  energies.  The  investigation  of  the  angular  dependence  of  the  cross 
section  showed  that  the  largest  cross  sections  occured  at  0°  in  the 
laboratory.  Therefore,  we  decided  that  0°  would  be  used  in  the 
feasibility  study  of  boron  profiling  at  Ep^3.2  MeV.  Further  we 
decided  that  the  B  metal  target  would  be  used  as  a  standard  in  the 
profiling  experiments. 


TABLE  I.  Published  Cross  Section  for  the  B(p,n)  C  Reaction. 


VI.  PROFILING 


A.  Theory. 

The  purpose  of  this  chapter  is  to  discuss  the  feasibility  of 
using  the  ^B(p,n)^C  reaction  to  determine  concentrations  of  boron 
in  solids.  As  used  here,  profiling  means  determining  concentrations  of 
boron  as  a  function  of  depth  below  the  surface  of  a  target  of  other 
wise  known  composition.  The  profiling  is  accomplished  by  comparing 
neutron  time-of-f light  spectra  from  two  targets,  one  of  pure  boron  and 
the  other  of  an  unknown  concentration. 

One  example  of  neutron  TOF  spectra  for  a  profiling  experiment  is 
shown  in  Fig.  12.  The  two  TOF  spectra  were  obtained  with  the  same 
experimental  set  up  used  in  the  cross  section  determination.  The 
upper  TOF  spectrum  (squares)  was  obtained  from  a  pure  boron  metal  tar¬ 
get  with  an  incident  proton  energy  of  4.517  MeV,  flight  path  of  2.197  m 

and  an  integrated  beam  charge  of  200  uC.  The  lower  TOF  spectrum 

18  18 

(circles)  was  from  a  boro-silicate  glass  target  with  the  0(p,n)  F 
TOF  spectrum  from  pure  Si09  (Fig.  4)  subtracted  from  it.  The  lower 
spectrum  was  obtained  with  the  same  incident  proton  energy  and  flight 
path,  but  an  integrated  charge  of  100  uC.  Statistical  fluctuations  are 
more  noticable  in  the  lower  spectrum  because  the  pure  boron  spectrum 
has  been  divided  by  two  to  normalize  the  integrated  charge  to  100  uC- 
However,  it  is  very  evident  that  the  boro-silicate  glass  TOF  spectrum 
shows  a  strong  resemblance  to  the  pure  boron  TOF  spectrum.  The  two 
spectra  have  the  same  general  characteristics. 
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The  theory  and  application  of  profiling  elements  by  neutron  time- 

10  22  23 

of-f light  has  been  studied  in  the  past.  ’  The  method  used  in 

this  thesis  was  based  on  the  article  by  Overley  and  Lefevre.  The 
concentration  of  boron  atoms  as  a  function  of  depth  can  be  determined 
in  the  following  manner.  Let  be  the  atomic  stopping  cross  section 
for  an  atom  of  type  i.  Assuming  the  sample  of  unknown  concentration 
of  boron  consists  of  boron  in  a  host  material,  one  can  then  define  '4>„ 
as  the  molecular  stopping  cross  section  for  the  host  material, 
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where  n^  is  the  number  density  of  host  molecules  and  n  is  the  number 
density  of  atoms  of  type  i  in  the  host  material.  Eq.  6.1  is  related 
to  Eq.  2.3  where  the  stopping  cross  section  per  atom  is, 
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where  n^  is  the  number  density  of  11B  atoms  and  n1Q  is  the  number 
density  of  10B  atoms.  Expressing  the  two  TOF  spectra,  taken  at  the 
same  laboratory  angle,  in  terms  of  their  differential  cross  section 
(Eq.  2.2)  one  obtains, 
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where  the  a'  indicate  the  number  densities  of  boron  atoms  in  the  host 
material.  Taking  the  ratio  of  the  two  TOF  spectra  as  expressed  above 
and  assuming  experimental  parameters  are  identical,  one  obtains  the 
atomic  concentration  of  boron  as  a  function  of  depth  x. 


nB(x) _ 

nB(x)  +  n^ (x) 


^Ndet(tV1  +  x’> 
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where  n^Cx)  is  the  number  density  of  all  boron  atoms  and  x  =  ^o^ll' 

For  a  pure  boron  sample  x  is  equal  to  0.2/0. 8  and  for  an  implanted 

sample  of  ^B  atoms  x '  is  equal  to  2ero. 

Integrating  the  stopping  cross  section  for  the  implanted  sample 

22 

one  can  calculate  the  depth  x. 
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where  t  in  is  the  time  corresponding  to  the  flight  time  of  neutrons 

produced  at  the  front  surface  of  the  target,  dE^/dE^  is  obtained  as 

discussed  in  Chapter  IV  and  dEn/dt  is  obtained  from  the  neutron  flight 

time  and  energy  relationships  mentioned  briefly  in  Chapter  II. 

The  depth  and  corresponding  concentration  (Eq.  6.5  and  6.6)  were 

evaluated  by  numerical  approximations.  The  technique  outlined  by 
22 

Overley  and  Lefevre  was  used  in  this  thesis.  Basically,  this  involved 
taking  a  depth  x  and  a  depth  interval  Ax  and  calculating  the  ratio  of 
the  number  of  counts  in  the  corresponding  time  interval  in  the  TOF  spectra. 
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B.  Resolution  Analysis. 


The  resolution  as  defined  by  Overley  et  al.^  is 

Ax  in  the  depth  at  which  the  reaction  occurs. 

Ax  -  AE  (En  i|>  )  L  , 

P  11 


the  uncertainty 
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where  AE^  is  a  proton  energy  interval  and  n ^  is  as  defined  previously. 

Assuming  the  front  surface  of  the  target  to  be  smooth,  there  are  three 

main  effects  that  can  contribute  to  the  resolution:  incident  beam 

energy  spread,  energy  straggling,  and  finite  time  resolution  which 

gives  rise  to  a  spread  in  neutron  energy  and  hence  in  proton  energy. 

Klystron  beam  bunching  introduces  a  5  KeV  spread  in  incident  beam 

energy.  For  a  S10,,  target  with  Ep  "  4.5  MeV  and  a  Si  target  with 

Ep  -  3.2  MeV  the  associated  depth  uncertainties  are  0.26  urn  and  0.31  urn 

respectively.  The  values  of  En^i^  in  Eq.  6.7  were  calculated  with 

13 

Eq.  2.4  in  units  of  MeV/um  from  Janni's  values.  For  a  SiO,,  target 

2"’  -3 

with  a  molecular  number  density  of  2.326  x  10  “cm  the  calculation 

-2  2 

constants  of  Eq.  2.4  were  -  1.776  x  10  MeV“/um  and  C.,  ■  2.568. 

22 

The  values  for  a  Si  target  with  a  number  density  of  4.997  x  10 
atomis  cm  ^  were,  *  1.738  x  10  “  MeV“/um  and  C.,  ■  2.322. 

Straggling  has  an  effect  only  below  the  surface  of  the  target.  As 

the  proton  penetrates  further  into  the  target  the  straggling  becomes 

"*4 

more  important.  Utilizing  Bohr's  expression  for  energy  straggling"  one 
can  solve  for  the  uncertainty  Ax^  (FWHM)  due  to  straggling, 
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where  Ax,  and  x  are  in  microns  and  n^^  is  evaluated  at  the  energy 

the  proton  has  at  depth  x.  The  constant  c  is  equal  to 
4 

2.355  [4ne  2(z^n^)]‘  where  z ^  is  the  atomic  number  of  atom  i. 

Differentiation  of  the  basic  relationship  between  neutron  energy 
3nd  flight  time  leads  to  the  depth  uncertainty  Ax^  (FWHM)  due  to  time 
resolution  At, 
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where  X  is  the  flight  path,  m^  is  the  mass  of  the  neutron  and  At  is 
the  time  resolution  (M..5  nsec).  The  stopping  power  En^'^  is  again 
evaluated  at  the  energy  the  proton  has  when  it  produces  a  neutron  of 
energy  E^.  Eq.  6.9  assumes  AE^  *  AE^  as  discussed  in  Chapter  V.  Unlike 
the  other  two  effects,  Axt  can  be  minimized  experimentally  by  increas¬ 
ing  X  or  decreasing  the  neutron  energy.  Neutron  energy  can  be  decreased 
by  decreasing  the  incident  proton  energy,  so  long  as  the  proton  energy 
at  any  given  depth  remains  above  threshold. 

The  results  of  quadrat ically  combining  the  three  contributions  to 
depth  resolution  for  a  SiO,,  and  Si  target  are  shown  in  Fig.  13.  The 
incident  proton  energies  and  flight  paths  are  indicated  below  the  respec¬ 
tive  curves.  The  curves  end  at  that  depth  where  the  proton  energy  has 
reached  threshold,  10.1  um  for  Si  and  79.6  um  for  SiO,,.  When  the 
curve  rises  with  depth  straggling  is  dominant  and  where  the  curve  falls 
with  depth  flight  time  resolution  dominates. 
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A  computer  code  for  calculating  lithium  concentrations  in 
12 

niobium  was  moditied  for  boron  in  silicon  and  for  boron  in  silicon 

dioxide.  (The  program  for  silicon,  SIBPRP  is  in  Appendix  F.)  The 
second  code  was  applied  to  the  two  spectra  of  Fig.  12  to  calculate  the 
concentration  of  boron  atoms  as  a  function  of  depth  in  the  boro- 
silicate  glass  sample.  The  results  of  the  calculation  are  shown  in 
Fig.  14.  As  expected,  the  concentration  of  boron  is  constant  at  about 
5  atoms  per  100  SiO,  molecules.  The  error  bars  are  based  on  the 
statistical  uncertainty  in  the  number  of  counts  in  the  neutron  TOF 
spectra.  The  negative  concentration  values  are  due  to  statistical 
fluctuations  in  the  two  TOF  spectra.  A  front  surface  resolution  of 
•^4  ym  is  experimentally  obtained.  This  is  in  contrast  to  a  calculated 
resolution  of  2.6  um.  One  reason  the  experimental  resolution  for  the 
boro-silicate  target  might  be  worse  than  calculated  could  be  due  to 
Impurities  near  the  target  front  surface.  This  particular  target  had 
been  used  for  a  number  of  years  to  align  the  accelerator's  charged 
particle  beam,  as  such  has  been  bombarded  with  a  lot  of  protons  and  a 
particles  over  the  years. 

C.  Sensitivity  Analysis. 

To  determine  the  minimum  concentration  of  boron  atoms  one  can 
detect  using  the  ^B(p,n)^C  reaction,  the  sensitivity  of  the  method 
must  be  investigated.  The  sensitivity,  defined  as  the  number  of  boron 
atoms  per  unit  area  n  Ax,  which  can  be  detected  can  be  determined  from 
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Eq.  2.2  and  Eq.  6.7  as. 


Overley  et  al.  solved  for  the  limit  of  the  sensitivity  in  the 

following  manner.  If  the  areal  density  of  boron  atoms  is  to  be 

determined  to  ±10%,  then  the  number  of  counts  N,  (t)  At  must  be 

det 

greater  than  100.  If  one  assumes  a  beam  current  of  1.5  pA  is  incident 

on  the  target  for  1  hour,  that  a  detector  11.43  cm  in  diameter  with  an 

efficiency  of  30%  is  located  1.593  m  from  the  target,  and  that  the 

reaction  cross  section  is  9.5  mb/sr,  then  the  minimum  boron  content 

14  -2 

which  can  be  measured  is  n  Ax  *  2.8  x  10  cm  .  Thus,  if  it  is  desired 

D 

to  determine  the  boron  content  to  ±10%  in  a  0.1  pm  depth  interval,  the 
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density  of  boron  atoms  must  be  greater  than  2.3  x  10  cm  .  For  boron 
in  otherwise  pure  Si  this  represents  560  atomic  parts  per  million  (ppm) 
and  in  otherwise  pure  Si0?  this  represents  a  concentration  of  1200 
boron  atoms  per  million  SiO,  molecules. 

In  practice  one  cannot  expect  to  obtain  the  above  calculated  lower 
limit  of  sensitivity  because  of  several  types  of  background  effects. 
There  are  three  main  background  effects  that  can  reduce  the  sensitivity. 

The  first  is  a  time  uncorrelated  background  cosmic  rays,  natural  radio¬ 
activity,  and  from  roora-thermalized  neutrons  inducing  (n,y)  reactions. 
These  effects  can  be  minimized  by  building  a  better  shielding  system  and 
by  using  a  pulse  shape  discriminator  to  distinguish  between  neutron  and 
Y  ray  signals  from  the  scintillator.  Other  sources  of  background  are 
time-correlated.  These  can  result,  for  example,  from  (p,n)  reactions 
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from  light  nuclei  target  constituents  with  low  (p,n)  threshold  values. 

18 

An  example  of  this  is  the  boro-silicate  glass  target  where  0  was 
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present  and  neutrons  from  the  0(p,n)  F  reaction  had  to  be  subtracted. 
The  third  main  effect  is  from  scattered  neutrons.  If  boron  atoms  are 
localized  at  one  depth  in  the  target,  then  there  will  be  neutrons 
emitted  with  a  very  well  defined  energy.  These  neutrons  can  then  be 
scattered  in  the  room  and  be  detected  at  lower  energies,  thereby  pro¬ 
ducing  a  low  energy  tail  in  the  TOF  spectrum.  This  tail  then  tends  to 
mask  the  effect  of  any  other  boron  concentration  at  greater  depths, 
thereby  reducing  the  sensitivity. 

To  investigate  the  sensitivity  experimentally,  basically  the 
same  experimental  set  up  was  employed,  but  with  the  detection  system 
used  by  Lunnon.^  To  reduce  time-uncorrelated  background  effects  the 
scintillator  system  was  shielded  and  pulse  shape  discrimination  was 
used.  A  block  diagram  of  the  electronics  is  shown  in  Fig.  15.  The 
scintillator  was  a  11.43  cm  diameter  by  2.54  cm  thick  NE  218  (Nuclear 
Enterprise)  liquid  scintillator.  This  scintillator  was  chosen  by 
Lunnon  because  of  its  good  discrimination  between  neutrons  and  y  rays. 
The  scintillator  shielding  consisted  mainly  of  a  1  meter  long  by  1 
meter  diameter  doughnut  shaped  tank  filled  with  boric  acid  and  water. 

The  scintillator  was  mounted  inside  the  doughnut  hole  and  shielded  by 
lead  blocks.  The  pulse  shape  discriminator  used  in  this  investigation 
is  described  in  detail  by  Lunnon. ^  The  neutron  window  was  set  with 
a  single  channel  analyzer  (SCA)  using  an  Am-Be  neutron  source.  The 
SCA  was  adjusted  to  eliminate  as  many y  rays  as  possible  without  losing 


too  many  neutrons.  The  resulting  signal  was  placed  in  coincidence 

with  the  SCA  setting  a  low  and  a  high  pulse  height  bias.  The  low  bias 
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pulse  height  spectrum  was  set  to  1/5  of  the  Am  60-KeV  y  ray  full 

energy  peak.  The  timing  bias  was  adjusted  to  the  low  bias  pulse 
height  spectrum  as  described  in  Chapter  II.  The  high  bias  was  used  to 
reduce  background  effects  as  much  as  possible,  and  was  adjusted  each 
time  the  incident  proton  energy  was  changed.  The  high  bias  was  set 
to  cut  off  pulses  greater  than  those  produced  by  1.69  MeV  or  0.31  MeV 
neutrons  for  incident  proton  energies  of  4.517  MeV  or  3.216  MeV  respec¬ 
tively. 

In  order  to  maximize  the  sensitivity  limits  of  this  experiment 
one  must  have  a  high  cross  section  and  a  large  solid  angle  (small  flight 
path) .  To  maximize  resolution  one  should  have  low  neutron  energy  and 
long  flight  path.  Thus  to  give  a  reasonable  depth  resolution  and  the 
best  sensitivity  an  incident  proton  energy  of  3.2  MeV  was  used  and 
neutrons  were  observed  at  0°  with  a  1.593  m  flight  path.  This  gave  a 
maximum  neutron  energy  of  0.31  MeV,  a  fairly  high  cross  section  of 
9.5  mb/sr  and  a  minimal  flight  path  consistant  with  the  constraints 
imposed  by  the  shielding.  Neutron  TOF  spectra  were  taken  either  at  1 

or  1.5  PA  beam  current  for  around  one  hour  on  a  boron- implanted  <100> 

* 

silicon  target  doped  with  either  As  or  P.  Each  targets'  TOF  spectrum 
was  taken  for  a  total  of  5  mC  integrated  beam  charge.  The  targets  were 
silicon  wafers  implanted  at  200  KeV  with  1  x  10^  atoms  of  cm 

it 

The  targets  were  supplied  by  Dr.  Dean  Casey,  Manager,  Material 
Research  Laboratory,  Tektronic,  Inc. 
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(isotopically  separated).  An  example  of  a  neutron  TOF  spectrum  is 
shown  in  Fig.  16.  The  small  peak  ("^120  counts)  near  channel  420  cor¬ 
responds  to  neutrons  emitted  from  the  implanted  boron.  The  large 
Y  ray  peak  lies  beyond  channel  650.  Again  the  channel  width  was 
approximately  0.97  nsec /channel. 

Immediately  following  the  boron-implanted  target  run,  a  background 
spectrum  was  taken  under  the  same  experimental  conditions.  The  boron- 
implanted  Si  target  was  turned  over  and  the  back  of  the  target  was 
bombarded  with  protons.  Fig.  17  shows  the  background  TOF  spectrum 
obtained  for  a  total  integrated  charge  of  5  mC.  The  flat  background 
below  channel  550  indicates  that  the  peak  in  Fig.  16  near  channel  420 
was  indeed  due  to  the  implanted  boron.  The  peaks  in  both  Figs.  16  and 
17  are  basically  the  same  in  channels  beyond  550.  The  two  smaller 
peaks  in  Figs.  16  and  17  near  channel  600  are  thought  to  be  a  result  of 
higher  beam  mass  or  gating  difficulties  as  described  in  Chapter  III. 

The  peaks  in  channels  beyond  700  are  thought  to  be  due  to  the  24  MHz 
bunching  rate  for  the  proton  beam.  Careful  investigation  of  these 
peaks  reveals  that  the  peak  near  channel  710  is  indeed  separated  by 
40  nsec  from  the  prcmpt  >  rav.  Further,  one  can  see  in  Fig.  17  another 
peak  40  nsec  earlier  near  channel  750.  These  peaks  are  very  small  as 
compared  to  the  main  y  ray  peak  and  are  seen  only  in  runs  of  long  time 
duration. 

The  spectrum  of  Fig.  16  was  corrected  for  background  (Fig.  17)  and 
was  analyzed  in  much  the  same  way  as  the  boro-silicate  glass  data.  A 
TOF  spectrum  from  a  pure  B  metal  target,  obtained  under  the  same  experi- 
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mental  conditions  was  used  as  the  standard.  The  program  SIBPRP 
(Appendix  F)  was  used  to  obtain  the  concentration  of  boron  atoms 
as  a  function  of  depth.  Three  Si  targets  were  investigated,  two  As 
doped  and  one  P  doped.  Their  respective  concentration  profiles  are 
shown  in  Fig.  18.  The  data  points  at  less  than  zero  depth  are  due 
to  calculation  of  the  concentration  in  front  of  the  target  and  give  an 
indication  of  the  background  fluctuations.  If  the  backgrounds  were 
identical  the  concentration  would  be  zero  in  this  region.  The  solid 
curves  in  Fig.  18  are  drawn  only  to  aid  the  reader's  eye  in  deter¬ 
mining  the  general  shape  of  the  peak. 

In  general  all  three  targets  show  a  peak  near  a  depth  of  1  urn  and 
have  a  maximum  ^  B  concentration  of  ^225  ppm  in  a  0.1  ym  depth  interval. 
The  actual  concentration  is  undoubtedly  much  higher.  It  appears  to  be 

decreased  by  resolution  effects.  Integrating  the  concentration  in  each 

11  14 

of  the  three  peaks  produces  areal  B  densities  of  9.94  ±  1.32  x  10 

cm  9.25  ±  1.39  x  10^  cm  and  8.33  ±  1.27  x  10^  cm  for  the  top, 

middle  and  bottom  profiles  respectively.  The  error  represents  the 

combined  uncertainties  in  the  terms  of  Eq.  6.5  assuming  the  stopping 

cross  sections  have  an  uncertainty  of  ^6%.  The  combined  uncertainty 

was  ^15%  due  mainly  to  counting  statistics.  The  experimental 

areal  densities  for  the  two  As-doped  targets  agree  with  the  implanted 
15  -2 

densities  of  10  cm  to  within  the  experimental  uncertainty.  The 
density  for  the  P-doped  target  lies  just  outside  the  uncertainty 
limit.  If  the  sensitivity  limit  of  concentration  is  calculated  for  a 
measurement  with  15%  precision,  one  obtains  199  atomic  ppm.  This  is 
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Fig.  18.  Concentration  of  boron  atoms  versus  depth  for  silicon 


implanted  with  boron.  The  silicon  targets  were  doped 
with  either  As  or  P  as  indicated.  Isotopicallv  separ¬ 
ated  B  was  implanted  in  the  Si  at  an  energy  of  200 
KeV.  The  concentration  was  calculated  assuming  a 
natural  abundance  of  boron  isotopes.  The  ordinate  must 
therefore  be  multiplied  by  0.8  since  only  is  present. 
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close  Co  the  experimental  value  of  225  ppm.  An  analysis  of  the  depth 
resolution  of  Fig.  IS  shows  a  FWHM  of  0.8  inn.  This  is  in  contrast  to 
a  calculated  depth  resolution  of  0.4  um  (Fig.  13).  The  difference  in 
depth  resolution  could  be  an  indication  of  the  limit  of  the  technique, 
hut  more  likely  it  is  due  to  an  underestimation  of  the  straggling 
effects  as  calculated  in  Eq.  b.S.  Bohr's  t'ormuLa  is  a  fairly  crude 
approximation  of  straggling  due  only  to  electronic  stopping. 

D.  Conclusions. 

The  final  objective  of  this  thesis  has  been  met.  The  ^B(p,n)^C 
reaction  can  be  used  to  determine  boron  atomic  concentration  to  v225  ppm 
with  a  precision  of  157..  The  theoretical  limit  of  resolution  was  not 
achieved  in  either  of  the  two  experiments,  0.35  pm  for  Si.  and  2.6  pm 
for  S10,.  However,  this  could  be  due  to  target  front  surface  problems 
and  straggling  estimations.  We  find  that  200  KoV  ^B  atoms  implanted 
in  silicon  reside  at  a  depth  of  vl  um  with  a  FWHM  of  0.8  pm. 

Tl\e  ^B  atom  implantation  depth  in  Si  has  been  studied  in  the  past 

25  26  >7- 

using  secondary  ion  mass  spectroscopy ,  ’  differential  capacitance,” 

and  Hall  effect  and  sheet-resistivitv^’^  techniques.  Boron  implan¬ 
tation  at  three  silicon  crystalline  configurations  have  been 
Investigated;  •-110''  which  Involves  channeling  and  <111^  and  ■•100''  where 
the  effect  of  channeling  is  reduced.  Further  studies  have  been  made 
on  amorphous  silicon  where  there  is  no  crystalline  structure  present. 
Implanting  ^B  atoms  along  the  >•110''  configuration,  maximizes  the 


29 
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effects  of  channeling  and  a  depth  of  1.7  t  0.2  pm  with  a  FWHM  of  0.6 

pm  at  an  implantation  energy  of  200  KeV  was  reported  by  Lecrosnier  et 
27 

al.  These  authors  also  investigated  the  <111>  cont iguration  at  the 

28  29 

same  energy  and  obtained  a  depth  of  0.8  pm.  Davis  ’  and  Lecrosnier 
27 

et  al.  have  studied  the  mis-alignment  of  the  crystal  off  the  <111' 

axis  with  ^B  atom  implantation  energies  of  200  KeV  and  obtained  values 

of  0.66  pm  with  a  FWHM  of  0.6  pm  and  0.6  pm,  respectively.  Ohmura 

et  al.^  studied  the  mis-alignment  of  the  <100'  crystal  axis 

and  further  reduced  channeling  effects  by  growing  a  thin  film  of 

SiO.,  over  the  crystal  surface.  At  an  implantation  energy  of 

200  KeV,  atoms  were  reported  to  have  a  depth  of  0.58  pm  with  a 

FWHM  of  0.25  um.  However,  it  should  be  noted  that  the  samples  are 

28 

still  crystalline  and  effects  of  channeling  can  still  be  there. 

Amorphous  silicon  was  also  studied  and  values  of  0.5  ±  0.1  pm  have 
26  29  32 

been  reported."  ’  ’ 

We  could  not  locate  a  report  specifically  on  the  unobstructed 
Si  <100>  configuration  at  an  energy  of  200  KeV,  but  conclusions  based 
on  the  above  data  can  be  drawn.  One  would  expect  the  <100'  and  ^111' 
conf iguration  of  silicon  to  have  an  implantation  depth  of  atoms, 
at  the  same  energy,  to  be  greater  than  that  in  the  amorphous  sample 
and  less  than  that  observed  in  the  <100'  configuration.  Based  on  this 
conclusion,  which  holds  true  for  the  <111'  configuration,  our  results 
of  the  <100'  oriented  silicon  sample  seems  reasonably  consistent. 
Further,  Wittmaack  et  al."  have  reported  that  the  depth  distribution 
is  not  gaussian  in  shape  at  implantation  energies  greater  than  40  KeV. 
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This  also  fit  our  data,  as  can  be  seen  in  Fig.  18.  The  reason  for  Che 
non-gaussian  shape  is  being  investigated. 

The  ^B(p,n)^C  reaction  offers  a  non-destructive  technique  for 
determining  boron  atomic  concentration  as  a  function  of  depth.  How¬ 
ever,  this  technique  probably  is  too  limited  in  resolution  and 
sensitivity  to  be  of  practical  use  to  the  semiconductor  industry,  where 
depth  resolutions  of  less  than  0.1  urn  are  desirable. 

Throughout  this  profiling  investigation  the  main  limiting  factor 
was  in  the  counting  statistics  of  the  foreground  and  background  TOF 
spectra.  These  can  be  improved  by  counting  for  longer  times,  but  a 
factor  of  10  improvement  in  sensitivity  would  require  counting  for 
about  four  days,  and  the  resolution  would  still  be  limited  to  about 


0.35  um  FWHM. 
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OOOOOOOOOO  OOOOOOOO. 3  00000000000000 OOOO 

CM  C«I  Cj  "•  »C\  O  a-  'O  t*-  tt.  •>  O  *»i  CM  v  *i>  (»  »r*.  —  CM  <T  —  W>  O  —  M  V  U>  «I  N  I,  O  C  —  u  «T 

r*0«t— 'PitvooocO  p  —  —  v^O^jp.m'*’.  —  c  p  u.f‘“a.*'M?f^u  — ov— >  —  p  p 

-Omt-'f'^-.uOO  u  a  r-  vp  ^  »o  ^  Cv.  CM  Vw  iM  Cl  CM  CM  CV  C\,'  \  -  o  l  1#N  — 

r-  \&  ip  ip  wi  v*  vp  i*j  14,'  c\  lit  vf\  >n  «ii  if\  u>»  <cv  »rv  >r\  «in  ai  'it  >r»  CMr>  «v«\  m  «i  s 


o 


a 

U) 


v  ^  v)  Mi  P  o  —  C\  ^  v  yA  o  a  C  o  —  P1  «/  /i  i)r  a  o  O  —  C  -  «■  t’  r*  <v  u  o 
—  —  —  —  —  —  C.  CM  CM  CN.  CM  CM  C  CM  C  C-  »»  •»»  «•>  »  ^  V  M  U  «J  V  u  «l  u  v  v  >P 
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APPENDIX  C 


Program  TOF-NE 


c  tpf-n'-  src  time  nr  m-ICvt  to  neutron  enercv  T<’ANFrr°YATirN 
c  ppnG'»A'«j  < a*>  ictn )/(<■'.?. «mev*f°)  vf  fne°cv 

C  FOR  2  1/2"  FOINT.  (rm-FfiA) 

C  FV?E'>I*ENTAL  PARA"ETFNF  1  tof  C»TA  ON  TAPE  UNIT  *1,  vit«  EFFIC  IFNCv 
C  TA?LE  (IF  USED  IN  PLCCF  1 

C  NFUTPOV  ENERCY  EPEOTRUr  OUTPUT  ON  TAPE  UNIT  42 
C 

C  NODI  FI  FT  7/11/7!  FOP  102A  CVANNEL  INPUT  FPEOTcUr  --  ON  C  A«  *<!  -  3 
C 

DCU2LF  “PFCISIQN  PTL.PTU 

DIMENSION  FFTAPLC  50) ,  NTOF(  1024), NNE(  102A)  ,T0F(4) 

00*VON  EFTAPL,  NTCF,T‘E 

data  “li,ul2,pl:/;u  tof-,?hne  7/ , i 1/71/ 

DATA  A,P,r,r/5HTvPF  ,5V1  c0e,5M  FFF:,5P  1.0,/ 

DATA  E,F,C/5H  2.0  .S^FCP  T,5HA?LF  / 

TATA  n, P, S, T, U/SHTEN  :,5HINPUT,5M  TOF  ,5WFV#  t,5HP0N#:/ 

DATA  P I , ?2, 53/5H  ENNA.5HX  :  .5WDCNE  / 

DATA  P!,E0/3.I4I6,«3«.5G5/ 

C 

r 

C  IDENTIFY  PRC CP  AY 

CALL  A*CI1(“L1,3) 

C 

r  A ?K  EFFICIENCV  to  ?E  UFEDt  FFFr  1.0  TYPE  1  .(it  =  rIT( EFTAPL)  Tn°E  " 
r  EFTAPL  (EFFICIF^fv  TAFLE  IF  F.°.  V/  FNLOW  1  ?t  FNT°'' ,  ENI'lC  END 
100  CALL  A  EOI I ( A , 7 ) 

CALL  ANVI NP( EFFFLG) 

*  :  EFFcLO 
50  TO  (200,300),  K 

ro  to  ioo 
c  Frcrrrr'<rv  :  i  .0 

200  FFF  :  1  .0 

CO  TO  400 

C  EFFICIENCY  :  EFTADL  ON  PUOr»  1,  PFAT  »ND  FET  POCCPA"  FLA G 
300  CALL  DTAPE(  1 ,  I  ,0,256 , EFTA °L(  D) 

C  NEUTRON  ENERGY  INCRFVENT  IN  *£'/ 

400  CALL  AFC! I ( ",  I ) 

CALL  ANviMO(r>F) 

C  DATA  TA^F  f’JIT'V  •  (2  TC  IS)  AND  1024  CHNL  .  °EC.IOM  4(1  TC  * ) 

500  CALL  ASCIKP.J) 

CALL  ANYIVP(FNI) 

600  CALL  AfCIKl',1) 

CALL  ANYI V°(PGN) 

J3  :  20  .•  (  CWI  *  1  ,H4,»(“GN-2.) 

C  TF’O  TOF  l  NE  eEGI ON? 

DO  POO  1:1 ,1024 
NTOF(  I)  :  0 
POO  NNF ( I )  :  0 

C  rAYI'”.'N  NEUTRON  ENERGY  IN  »EV 
CALL  AFCIK2I.2) 

CALL  ANN  I NR( FN  AX  > 

J  :  FNAt/”'E 

"  LOAD  Tu E  TCr  FBFOTPU" 

"ALL  DTAPF( 1,JP, 0,1074, NTOF< 1)) 

DO  not,  l;|,4 

I )  :  \'TOF(!> 

t/j 

g 


r*c* 

x  j  uu  a 


70 


-500  TCF(I)  :  Miff  I ) 

C  f|vPFn  I  TF'lTAL  °APAFETEPSt  rT<  FSE0/OMAVNEL)  ,  GATA(  '•FA'JVEl  »), 

r  El( FLIGHT  PATH  I V  -JETERS),  tTOl'LC ORAPGF  t;i  *!  05O00l,L0,'D  ?>  , 

TT  r  T0E< 1 )/ 10000, 

GAf-*A  :  TO  F(  2)  / 1  00  . 

EL  :  TO F(  2  )  / 1 000 . 

UOOL'L  :  TOF(4)/IOO. 

C O'* EGA  r  PI  *(0.02  H/ELW  0 ,0.M  VEL) 

TE^O  :  GAFKA  *  3.24‘EL/TT 
F»U  :  D  F*  0 . 5 

C  GENERATE  THE  \T UTRON  ENERGY  SPErTPL'T 
ro  opr  isi,.! 

IPl  :  1  +  1 
XI  5  I 
EN  :  TE*»I 

C  EFFICIFNC''  OF  PETEOTOR 

GO  TO  (020,510),  x 

C  'oiO'  IF  EFrIC!INCY  FROF  TAPLE  (Ol'BIC  rIT) 
olO  CALL  rI  T?  (  E'J,  EFF) 

IF(EFF.LE.O.O)  GO  TO  Poo 
o20  FML  :  E'W 

FNl-  :  OF* (*I+0.5) 

C  RELATIVISTIC  CALCULATION  OF  HISTrf.PA*  FOUNT A’lES 
GL  =  ENL/EO 
Git  r  FNU/EO 

PTL  :  2  .*  GL*  (  l  ,-3  ,/J2  ,»0L*  ( I  .  *  A  .  /  3  ,*  H. ) ) 

PTU  i  2,*GU*(l.-3./2  .•  "1‘*  (  1  .  *A  ,/G  ,*otj)  ) 

PETAL  :  "SooX( ? TL) 

PFTAU  :  DFnPT( PTU) 

TL  :  EL/(  .25<3  752*PFTAL) 

Tl!  i  EL/  (  ,?o<>7PJ*oETAL') 

OIL  r  TEup-TL/TT 
OIL'  z  7ER0-  TtJ/PT 
C  CALCULATE  THE  HISTOGRAM 

11  t  OIL*. 5 

12  :  CIl'+,5 
DELTA  :  OR’- OIL 

r  z  12-n 

Sl'Y  :  0.0 
IF(I1)555, 505,1 

1  IF(0ELTA-1 ,)2,2,3 

2  IF(F)A,4,2 

2  I  EC  F- 1)4,  7, 5 

<  Fi  :  t'ft.T* 

T?  z  0.0 

GO  TC  e 

5  N1  i  1 1  + 1 
F2  r  12-1 

00  G  Lr-M  ,  V2 

ICH  :  L+l 

xrATA  :  NTCF(ICH) 

6  SUV  :  FUv<-vDATA 

7  »It  :  II 
*1?  r  12 

FI  r  0 , 5-  (  OH  -x I  1 ) 

F2  :  0.5+OIU-XJ2 
®  I  ON  :  1 1  ♦  I 


Si 
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TCATA  :  VTOFUCM) 

CUT  :  F|  » "CATA 

1CH  :  12*1 

*CATA  :  NTCF(ICV) 

SI"*  :  FlT*F?»rPATA 

C  CALCL'LA  TICN  OF  «N*  1 0»<s  /  <  *  INC  .  PAPTICLFSaCEV*?? 1 

NNf(tPl)  -  SLY*  l .602E-A/(UC0LT*?F*CCrrGA*FFF) 

CALL  LITESC!) 
oor  CONTINUE 

C  VIP  ITT  OUT  NF'JTPCS*  ENERGY  SPECTPUT  ON  UNIT  »? 

CALL  ASCII (T, 1 1 
CALL  AimVPtFVJ) 

CALL  ASCII  CL',  I) 

CALL  ANY  I  V'P(  PGN) 

JP  :  20 ,«  <  ?'»'!  - 1 . 1 *4 .• (PGN-  2  . ) 

CALL  CTaPE<2,J3,I  ,1024,NNE<  IP 
C  TYPE  TONE* 

CALL  «SCI ! <? J, 1 > 

C  TYPING  1.  PESTAPTS,  2.  «ESTAPTS  I TM  EFF  AND  TEN  SET,  CTUE<>  STOPS 
CALL  ANYI'-Pm 
IT  :  y 

CO  TO  <100,5001,  IX 

STOP 

F\T 


JO  ooooonoon 
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APPENDIX  D 

Program  B11PNX,  Subroutines  KINSUB  and  FIT2 


91  1PNX  SBC 

FOR  1 1 B<  ° , M> 11C  RXN  ONLY,  CALC  TOP  TO  X-SECT 

CODIFIED  JAV  7?  Fro  OOBPN  P*N  *F0N  LI7AM  sd(>  by  SEaLCCY 

WHICH  '-'A?  ‘-on  Fim  rcc:*  LITPVl  «'ir  TFSTAT  -  OF?  TUEFIS. 

CFl.ri'LPTFS  SIGNA  IV  NICRCFARNF/S'3 

STA  TI STIC?  0UT°L'T  AS  10 ya /SORTC CNTS 5  (**10»0 

FOP  2  1/2“  SCINTILLATOR 

FOP  VAPIAOLS  STOPPING  CROSS  SECTION  -  0202, PN,P  TGTS . 


C 


COUPLF  P°F0!?I0N  STL. -’2U,  ARGLOP.DSCRT 

DIMENSION  EFTA3L< 1024) ,  KTC  c( 102 A) .^SlGl 1020 , TPF( ° 1 

coupon  eftabl. ytof, "•sis 

COP-OM  /YINCON/Al, A2,A3,EPY,  "’P  ,TC,  THETA  ,  FNC  ,  FN*  ,  D,>’0 

fO^'-ON  /SCTCCN/EN,  SP 

DATA  A3.AC/3H  PE  .SHIMFV)/ 

TATA  A  A  ,  3S/3HSVA  :,5MRC-N*:/ 

PAT*  CC.PO.EE,  PP/?HINPUT,5H  PT1  ,5HOUTpU,5HT  PT2/ 

TATA  PI ,ONE  *Nll/3 . I  a ! 6,Pi 1 .A? I / 

TATA  XI ,y2,r3,vA,-'P/l .OOTS25, 1 1 .CC«330,1 .00?f?5,l 1 .011402,1./ 
DATA  ?l,?2,r/!.,S.,»2.760/ 

TATA  PCI  ,rP2,DP3/5MI  VPl'T,  5H  TGT  ,5HA  / 

TATA  PD4,P’'5,PP6,2D7/3Hlr!>2C,  5“*  2:,”9"P  3:,5W?N  4/ 

PATA  Prp,P0P,rCl(j/5H:3N-O,5H  5sI,5«VPL'T  / 

PATA  rPl 1 ,PP12,DP13 ,DC14/9HINPUT,9H  STOP , 5VPI "0  .3WCRCFS/ 

CATA  rCI5,PDI6,rCI 7/5H  SECT, 5HICM  (,5«A,o)  / 

rn  z  i 
a  i  :  x  l 
A2  z  X2 
A3:  ><3 

PO  r  o 

n N  :  ri*OMFArO 
TIP  :  XI V 
Eli  :  x3»0NEA  Ml1 
xET  :  0.31 1004 

TYOTO  :  rp*'Yp*Y2*4.*PI*53r.l 6F- 40 /( YFY* 2.5S67E-I2) 

CONST  :  t.S0?.E-l‘»«l.E+06/(l.E-06*I.F-24) 

CALL  PTAPFd,  1  ,0, 1024, EFTA=L<  11) 

TPANSFCPHATION  CHANNEL  WIDTH  (NFV) 

CALL  AfCII(Aa.P) 

CALL  A  NY ! vp( 0  E 1 

INPUT  STOPPING  CROSS  SECTION  VALUE? 

CALL  ASCI  I < PD  1,31 
CALL  ASCIKCC4.7) 

CALL  AVvinR(CCI) 


C2  :  0. 

IFICCI .EP.I .)  Cl : 3, 325F-2 1 
IFfCCl .EP.l .)  C?:2.P3« 

I  F(  CC  l  ,E“  .2 .1  Cl:1.644r-?1 
IFfCCI  .EC. 2.)  r2:?.SF? 

1  F(  CC1  . EC  .3 . 1  Cl  :3  ,r  12E-2 1 
IFfCCI  .Er.3.1  C2:2.?7S.o 
I  Ff'TI .E-.A.)  Cl : 4 . 7°3F*2 I 
If<CCl .fr.4.)  r2iC.B46 
I?< "Cl .Sr, 5. )  GO  TO  4 
00  TO  7 

CALL  ASCtKCPl  1,7) 


no  oooom  oo; 
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CAU  ¥NYl»IP(Ct,C2,Pl,r2,r;> 

C  INPUT  CECTAPE  SWITCH  A,  PEGION  • 

IF(C2.E*.0.>G0  TO  5 
7  CALL  A  SCI  I (f C, 2) 

CALL  A  PC 1 1 ( A  A  f I) 

CALL  ANYlNP(fVM) 

CALL  ASCI  I  (S  B« I ) 

CALL  ANYIMP(RGK) 

C  OUTPUT  CECTAPE  SWITCH  #,  PEGICN  A 

CALL  ASCI  I (EE, 2) 

CALL  ASCI I( A  A, 1  ) 

CALL  ANV I NP( SVO) 

CALL  ASCI  I (EE, 1  ) 

CALL  ANYINP(RGO) 

JBV  i  20,*(SWN-1.>  ♦  4.*<RGN-2.) 

J30  :  20.*(SV0-t . )+4 .* (PGO-2 . ) 

C  CALCULATE  PROTO “  THRE5M0LC  ENFRGY 

EPPIN  t  CXI+x2)*(-?>/X2 
KT¥IN  :  EPMIN/DE 

C  EERO  DATA  A NO  SIGPA  PAT^ICES 

CO  222  P: I, 1024 
«TOF<  !■')  ;  0 

XSI3(M>  i  0 

2  CONTINUE 

REAP  TOF  SPEPTP'JP 

WPITF  TRANSFCR^EC  SPECTRUP  ONTO  TECTAPC  UNIT  2. 

CALL  DTAPECt, JEN, 0, 1024, KTOF<m 
DO  224  IX  :  1,7 
TOr< IT>  :  XTOF(IX) 

4  CONTINUE 

°EAC  »AnAvATEPE  FROM  DATA  TAPE: DT< FLI GHT  TINE, »<ANOSEC) 
GANVAC CHANNEL  *> , EL( NEUTRON  FLIGHT  dATH,  PE TEpS5 , 

UCOULC INTEPGPATEO  PEA*!  CURRENT , VICRC  CPUL) ,TUETA( CTG) , 
Enax  < I  NCI  TENT  PPOTCN  ENEPGY,»EV) 

DT  :  TOF( I) /1 0000. 

GAP* A  i  TOF(2)/IOO. 

EL  i  TOF( 3 )/ 1000 . 

UCOUL  =  T0F(4)/|0. 

THETA  :  TOF( 5) / 1 00, 

ETA X  i  TOFC 6) / 1000 . 

CCP  :  TCF<  7) / 100. 

IF(CCR.EF.O.O)  CCP  i  1. 000 
XT¥AX  :  EPAX/EE 

TOF  SCATTERING  COPRECTI ON— SLOPI NG  LINE  HVHP  50NEFC/VETER 
PAX  CHANNEL  AND  PA X  COUNTS  INPUT  WITH  TOF  CATA 
PAXCHN  :  KTOFO?) 

I F(  PAXCHN.  EO  .0 )  GO  TO  230 
C NTV A v  :  XTOF(P> 

PINCMN  r  IOO.*EL/DT 
H1NCHV  s  ¥AXCVN  -  MINCHN 
IF(PI''CHN*  1 3.LE.0)  PINCWN  t  15 
CO  227  ICMNL  :  *1  f.'CHN, PA *CHN 
T  i  ^AXCVN  -  tCHNL 
T  :  T-CT/EL 

ICNT  s  (IOO.-T)  •  CNTNAX/IOO. 

IF(ICNT.LE.O)  GO  TO  227 
YTOF(ICVNL)  i  *TOF( ICHND  -  ICNT 


ooo  o  no 


K 
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227 

2J0 


CONTINUE 

D0*EGA  i  3. IAI6*<.G315/EL)*C.  0315/EL) 

JE^O  :  GAMMA ■*-3.34* EL/DT 
EPIC  :  EPPIM  -  PE/2. 

CALL  XINSU? 

ENU  s  ENK 

CALCULATE  TRANSFORM 
DC  POO  ''T:XT’*IN,XTMAX 
XT!  :  KT  »  I 

CALCULATE  the  energies 

XKsXT 
EPtPEM* 

EPU:DE*CTX*.5> 

ENLsFNU 
EPX  s  EP 
CALL  « INSUP 
EM  s  F.NX 
EPX  i  EPU 
CALL  X INSUP 
EMU  s  FNX 
RELATIVISTIC 
GI.  i  ENL/EO 
GU  s  ENU/EO 

^  PTL  s  2,*GL*(1 .-3. /2,*GL*< I .-A./3. *CL> ) 

B2U  t  2  ,*Gl>  (1,-3  ,/2.*Gl>  (1  , -a ,/J  ,*GU)  ) 

BETAL-PP^PTCPTL) 

9ETAl’rrscRTCF2U> 

TL:SL/(  .2«07S3*'»ETAL) 

TU:EL/<  .2P«7«»3»3ETAU) 

CILi-’E'O-TL/PT 

CIUs?ERO-TU/PT 

GET  SCATTERING  CORRECTION  FOR  GIVEN  EN,T**ETA  -CISPFGA»D 
FOP  3 1 1  °M  PXN. 

SR  :  |  . 

EFFICIENCY  NOV  defined 
CAU.  rIT2( EN,  FrF) 
ircEFF.LE.0.0)  GO  TO  800 

CALCULATE  PFA/rr  (NEV  C*  CM)  FCP  Pll  PRCX  °M  TGT  CHPP  U^ING 
EMPIRICAL  EXPRESSION  CALCULATED  FROM  .’ANN I  VALUES  FOR 
PO°ON  AND  NITROGEN. 

DEADX  :  Cl /EP  •  <  ALOGC  EP)  *  T2) 

HISTOGRAM 


CALCULATION  OF  HISTOGRA*'  BOUNDARIES 


- 

1  : 


y 

u 


THE 


CALCULATE 
IlsOIL*.8 
I2sCll'*,5 
DELTArCIU-CIL 
Pi  12  - 1 1 
SU*tO. 

IF< II >800,800,28 

I F<  DELTA- t . ) 1 ,  1 ,2 

IF(M)2,3,2 

IF(r- I  12,26,  6 

FI  :  DELTA 

F2s.(j 

GO  TO  27 

K:  11*1 

»'s  It- 1 

CO  200  Lsv, N 


1 

. 
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I  CM  r  L+! 

XDATA  r  XTOE(ICW) 

SUM  :  SUM  +  XDATA 
200  COMTI MUE 

26  XI  1  :  It 
XI2  :  12 

H  :  0.5  -  (C1L-XI1) 

E2  i  0.5  +  CIU  -  XI 2 

27  ICH  :  1 1 1 

XrATA  :  XTOE(ICH) 

SUM  i  SUP*  *  Et*XDATA 

ICH  :  12+ 1 

XDATA  I  *TOE(ICH) 

SUf!  :  SUM  +  E2*»rATA 
c  UNCERTAINTIES 

C  SOUP CE  SCATTERING  UNCERTAI NT1FS--  X  OF  INITIAL  C«TA  CORRECTION 

C  USE  CHANNEL  AS  AVERAGE  OF  TWO  LIMITS 

SSPR  i  0. 

ICHSS  :  (OIL  +  CIU)/2. 

TC  :  MAXCHN  -  ICHSS 
TC  :  TC*DT/EL 
IFCTC.LE.O.O)  GO  TO  735 
IF<TC.GF. 100.0)  GO  TO  733 
TC  r  (100.-TC)*CNTMAX/I00. 

SSPR  :  XTOF( ICHSS) 

SSPR  i  TC/(TO+SSPR) 

SSPR  •  $SPR*SSPR 

C  BORON  SCATTERING  RESONANCE  UNCERTAINTIES  AS  X 

733  XLRPR  :  A3S( 1 .00  -  SR) 

XLRPR  :  XLRPR+XLRPR 

C  COUNTING  STATISTICS  UNCERTAI NTV  AS  X 

STATPR  :  1.0/  SUM 

C  COMBINE  ANr  OUTPUT  UNCERTAINTIES  AS  X  •  IOTA 
XTtP  i  XT  1  +  400 

BTL  t  STATPR  ♦  0.3* (SSPR  +  XL®??) 

KSIG(XTIP)  :  t .CE+04*DSGRT<BTL) 

C  CALCULATE  SIG^A  AND  DEPOSIT  IN  XSIG  MATRIX 

FSIG(NTl)  :  CEADX*OOP*SR* SUM* CONST/ <UCOUL*EEE* DOME GA*CE) 

BOO  CONTINUE 

DO  850  XL: 1,6 
KSIG(XL)  :  XTOE(XL) 
e50  CONTINUE 

C  WRITE  SIGMA  ONTO  TECTAPE  UNIT  *2 

CALL  DTAPE<2, JBO, t, 1C24,KSIG< I)) 

GO  TC  5 
END 


PV 


S' 
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FI "SUP  C3£ 

Kt'*EvATTC?  PUS POUT! VE 

( VI TV  "AVEO  CCf’-'CV  *  "INCO"'  rnp  °A  ?AVF  Tf 5  PASSING) 

"CPI  FI  EC  21  PERT.  le  70 ,  CA? 

pvn.  pevimptov,  ptastort  ir!vr=rrrr  s/u/f6 
rpfl'icff'nrp  energies,  angles,  avc  c°css  sections 
3ETvEF*l  or  AT  LA?  SYFTFPS 
FEE  NCCLEA®  TATA  TArLFF,  CA  "ES  161, IT  rf'c  “10  TATI  CM 
Al,  A2,  A3  Apt  MA PEES  <  IV  an’)  OF  ITirSVT,  TARGET, 

AN'C  FWI  TTF2  PARTICLES 

FI  :  IMCITFMT  E"FRGV  Cl"  YEV)  IV  L«?  SYFTFY 
o  :  "  CF  OEAPTICN  (IV  REV) 

THETACT  AT  T"ETala?  A<»F  EY I  EPIC"  ANGLE  CIV  SEG^IES) 

E3r"  AT  f:la°  ape  ENERGIES  OF  F-ITTFO  PA°TIOLE  Cl*)  »F V) 
pati C  r  C*  C®OSS  SECTION  /  LAP  OPCFF  SFCTIO" 

"CPE  -  0:  INPUT  :  Al,  A2,  A3,  FI,  ",  T^ETAC" 

NODE  s  Os  INPUT  :  Al,  A2,  A3,  ?|  ,  T^FTALA?,  E3LA  ? 

"CPE  »  f-S  INPUT  :  a  A 2,  «3,  FI,  0,  TVFTALA? 

SUBROUTINE  vtm?uB 

CCUBLE  °R  EC  I S I C  v"  AO,  FCEL,  FLEC , FB , RT, AAEf 

por*?*!  /nimpp:vai,»2,  ab.fi  .’'.thetac.thftal, r?r», fflaf, 

1  PATIO, "OrF 
II  IF  CTCE)  13.  I’,  12 

1?  SIMLA?  i  FIN(T«ETal  «  0.017/1532*25) 

CO  PL  A  *  I  CCPC  TWFT»L  *  0  .0  1  74  532*25) 

13  AA"  r  A)  *  A2  -  A3 

IF  (TCF)  31,  21,  31 
21  P  :  (J. 

rc  25  I  :  I,  3 

AA  :  ajo  -  "/o?  | .<7? 

AAEF  :  Ai.aj.ri.rj^flC 

25  '  :  C ( A3*Aa )«E?LA3  *  (A|-A4)»FI  -  2. * CSC RTC  a a  EE) 

I  «0<3fLA3)  /  aj 

31  ET  r  El  +  2 

IF  (ET)  32,  32,  Al 

32  F3CT  :  *0. 

E2LA?  :  -Ci  ,000 
PATIO  :  -0. 

IF  ("orS.O.T.O)  THETAO  :  0.0 
PETU°N 

41  A4  :  A4~  -  "/P31  ,A7? 

CE'IC”  :  CAI*A2)  *  CA3-M4) 

»  :  II  •  Al* A4/I)Ef!PP 
?  :  a  *  a  3/ a  a 

C  r  (FT  *  "* At / A2 )  «  A2*s3/rF"0" 

IF  (C. t.f.0,0)  GO  TO  32 
E3  rr  t  C  *  A4/A3 
AC  r  A*0 

POC74C  :  PPOPTCAO) 

IF  ("PBS)  51,  61,  61 

51  f  I  VP"  :  ®IV(  TWfTar  .  0 .0  I  74 53 2*2 5) 

E3LAO;  T?cr  ♦  2.0* SO RTA C *CCS ( TM ETA 0*0 ,01 74532*25 ) 

EC EL  r  r3Pv/F3LA? 

ep'Laa  ;  eivC"  •  opOOTCFCEl) 

SI  PT  t  F3CT/n  -  *!vl  A?*f  JVLA" 

I F  C’T)  32,  62,  62 
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l 


r> 


r 


62  ROOT  :  CSCPTCPT ) 

If  (fCCE)  *>1,  05,  04 
e|  rC'iTIVL'E 

I?  CFI'ILAR)  02,02,03 
92  CCNT I  »!VF 

THFTAL  :  THETAC 
CO  TO  100 
03  CONTINUE 

FP  :  E3 LAP/? 

COSLA  5  0  r?CRT(ES>  -  ROOT 

TV  FT  A  L  :  oO.  *  ATA  N(  -COSLA  P/SI  ''LAS  )  /  0.01  745320?* 
GO  TO  100 
04  CONTIN'UF 

E3LAE:3«  (COSLA  9*P00T  >  »(CCSLA  ?+i>OOT> 

05  CONTINUE 

If  (SINL«?>  o«,  96,  07 
06  CONTINUE 

THFTAC  :  T«ETAL 
GO  TC  100 
0  7  CONTINUE 

ELEC  1  F3LA3/F 3 Or 

SI'T*  :  ?!  ML  a?  *  CEC0T(FLFC) 

COSC*  z  (F3LAP  -  330'’  -  ?)  *  0.5/CCPTAC 
THETA 0  :  90.  *  ATANC  -CCSCr/SP'CH)  /  0.0174532025 
100  CONTINUE 

DATI9  :  ff'OTAC* ?0CT/F3L»? 

RETURN 
ENT  ' 


£ 


1: 


r  FIT2 

SU3°0l'TI \'E  FIT2CEIN.SIG) 

PI  YENS  ION  S< 250) 

CC^rCM  S 

el  4  ?<n 

FC  :  0(2) 

IF(EIN-( EL*EC) >20,21,21 

20  SIGiO. 

RETU°N 

21  ENTRY:  1  .-t-((EIN-rL) /EC) 

Kr  FVTPY+2 

C  FIRST  2  FNTPIES  *RE  FLO*  AMC  EIVC, 

»Y:K 

TO: S( w*  I  ) 

>M:?m 
N?:SO'*l  > 

YJ: ?C  v+2) 

*:EMTrv-vv+3. 

A: ( - YO+3 .* v I -3  .•v2+v 3) 
P:C6,«V0-I5.»VI*I2.*Y2-3,»YS) 

C: ( -  I  I .• vO*IR .•YI*o.«7  2*2.»v3) 
SIC:Y0+X*(r*lr*(P*r*A))/5. 

30  RETURN 

EMC 


t 


oooooooooaooooo 
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APPENDIX  E 

Program  LSQFT3,  Subroutines  CBFIT8,  LEGPOL  and  MINV/ 


LSOFT3  SRC 

A  LFAST  SQUARES  PR0G7A*  TUAT  FIT?  A  SFRIES 
OF  l.EGENPPE  POL'^onIALS  TC  POPON  1KP.N)  ANGOLA  P  Cl  STPIFL'TI  ONS 
INPUT:  FIG-'A  <THET«LAP)  IN  YICPOBA?NS/F»,  NO.  Or  POLYNOMIALS 
(NPOL  ,  FAY  r  6),  NO.  CF  ANGLES  (NANGLF  ,  “A*  r  ID 
PE  IF  ENEPGY  INTERVAL  OF  THE  INPUT  CROSS  SECTION  IN  PE V. 

INPUT  CROSS  SECTION  FUST  PE  I"  AN  INOPEASINO  SECl'EVOF  OF 
ANGLES  EY  SWITCH  ANC  REGIO"  NO.  UPPFR  A  LCVER  El  LEVELS,  °YN 
MASS  i  o  VALUES  ARE  VARIAPLF.  MOPIFIEr  FPCM  LSSFIT  AUG  7?. 
TOTAL  Y  SECTION  '  JC(LG)  )  IS  CIVICFC  ?Y  10, 

LEG  COFF  :  ,I0(LI)*I0  +  U  AVC  EPPOP  J0<L2)*I0 
USES  NODI  FI  EC  CTFIT«  AUG  7F . 


C 

40 

C 

C 

C 

C 

549 


COUPON  JD<  5632) ,  «J0(  3072) , ELSC( 20,4),A<6,2) 

COPNON  /PINCOP/  A|,A2,A3,FAK,O,TC,THETA,Fvr,ENV,P,>«0 
DATA  r i t 02, C3/5H INPUT, 5w  SV,e,5«EG  / 

DATA  BltP2,33,R4/5HSV,RF,5HG  F0R.5H  OL'TP.SHUT  / 

TATA  B5,36,B7/5HAPRAY,5H  STA»,5MTF  AT/ 

DATA  S«,3R,?I0,311/5H  -SVI ,  5HTCH  ,5H  -PEG,5UI0N  / 

DATA  Oa  ,  C5 ,  Of  /5W30?0N ,  5N1  1  (  PN, 5U  )  °YN/ 

DATA  C7,rp,C9,C!0,0l  1/5«INPUT,5«  0  IF,5w  NO  0,5HHANGt,5M  TO  */ 
DATA  C12,CI3/5H1,Y2,,5HM3,<5  / 

TATA  on ,rC2,rC3,Cr4,r5/5WNPCL  .5HIKAY  ,5Hf),  N,5«ANGLE,5H  f'A'V 
DAT*  D6.C7/5N  1 11,,5H  DF  / 

DATA  38, P«, DIO, PI  l.^/SH  TO  E,5HI  LIY.5HITF  ( ,  5"304 , 4 , 5H5 1  >  / 

DATA  P12A/5H  FI  / 

DATA  ri3,ri4,Pl5.PlS,ri7/5MC0S(T,5HC)  PA,5ur,  FR, 5HPCC ,  ,5HtSEC7 
DATA  PI?,,'|Q,r20/5M,  0AL.5H  LEG  .5HYSFCT/ 

PATA  014,0|5,Clf,C|7/5HNANGL,5HF  T00.5HLARGF ,*w  FO*  / 

DATA  Cl S,ri«tC20,C21 .C22/5NEI  PA.5HHGE  ,5VTTv  0 , 5HUTPUT , ?M  FES7 
PATA  Cr3,C24,C25/5M°EP7  ,5Hf  0  : ,  5H  NO  ]/ 

TATA  vi , Y2 , v  3 , v  4 / 1  .007«25,ll  .00*350 , 1  .0086 65 , -2 .763/ 

GIVFS  »ry  4  ALLOWS  FOR  CHANCE  CF  YASS  4  C  VALUE. 

CALL  »SCII(C4,3) 

CALL  A  SC II(C7, 7) 

CALL  MNYlNPCAAI, AA2,  AA3,0C(Ttri  ) 

ALLOWS  FOR  TTV  OUTPUT  OF  Y.$FCTICM. 

CALL  ASri!(C20t6) 

CALL  YNYINPCAPATA.PI ,r2,P3,r4) 

I*I°UT  THE  LOCATION  OF  THE  CIPST  aNCLF  CF  THE  SERIFS. 

C»LL  ASCI  I (C 1 ,3) 

CALL  "'NYINP(S'J,?EC,P3,ri,r?f 
INPUT  THE  STAPTING  LOCATION  POP  OUTPUT  CAT« . 


C*LL  ASCI  1(31,4) 

CALL  M,tv I ;;P(  SVO.PrCO.PI  ,P2,PJ  ) 
INPUT  *  or  DCLVN0~I»LS,  *  OF  ANGLES 
CALL  ASCII(PDI,7) 

CALL  MNYI NP( YPrL, VANGLF.DE.PI ,P? 

NPOL:  YO0I.*,5 

WANGLE: VANGLF+.  5 

MCi  I 

A  I :  v  I 

A2r»? 

A3  :»3 


IN  SERIFS  4  THE  FNFRGY  INTF°VAL 


^  j ^ 
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c 


541 

C 


c 

e 

c 

c 

c 

501 


C 

601 

602 

543 

C 

54  4 
201 
20? 


«'tT4 

CHANGE  OF  YA SS  4  G  VALUE  IF  REQUESTED. 

IF(AA1  .EO.O.)  GO  TO  541 
A 1  z  AA1 
A2  :  AA2 
A3  -  A  A3 

p  :  COO 

CALL  *FCI I ( C7, 4) 

I CUr  z  311 

ALLOW?  ?C?  CHANGE  OF  ENERGY  LIMITS. 

CALL  A  SCI  I  (  25 ,  5) 

CALL  MNYINPCXEIL.XEIU, Cl  ,22,23) 

I  El  L  :  304 
I  El  U  z  451 

IF  (XE1L.ES.0.)  GO  TC  501 
IE1L  z  yEl  L+  .5 
I  El  U  z  XE1U+.5 

CALCULATION  OF  THE  NATRIX  REFERENCE  A. 

IE1L  4  IE1U  z  ENERGY  CF  Ml  LOWER  AND  UPPER  VPT  v-SECTION. 
CONSTANT?  ICC1  TO  ICC10  ARE  FI XF2  FOR  TWE  GIVEN  ENERGY  REGION. 


ICC1  z  IE1U  -  IF1L 
ICC2  z  2  +  ICC1 
ICC3  z  ICC2  +  1  +  ICC1 
ICOA  ;  ICC 3  +  1  ♦  ICC1 

i ccs  i  irn  ♦  i 

ICC6  z  IE1L  -  ICC 3  -1 
I C07  r  IE1L  -  ICC2  -  1 
ICC®  z  I  El  L  -  2 
ICC®  z  I025-IE1L 
ICC10  z  ?04<5-I  El  L 

CALCULATION  of  NATPT v  2!  PENSION  ope  tip. 
ICC1I  z  512-NANGLE 
ICCt?  z  512 
IC4  z  TCC4-2 

IF  (ICC4.GT.512)  GO  TO  601 
GO  TO  602 

ICC11  z  1024  *  NANGLE 
ICC 1 2  I  1024 

IF  (ICCtl .GT.5632)  GO  TO  543 

GO  TO  544 

CALL  A  FCI I ( 0 1 4 , 6 ) 

GO  TO  545 


EERO  THE  *AT®!y  AND  RE»20NE  ANCLE  *T  «  TI*E  (/-SECTION, r5®''r)  , 


®A2  z  3. 1  4 1 5c /i po. 

re  201  I z  1  , 5  63  2 
«.'2(  I )  zO 

20  20?  Iz  1,3072 
J0( !)zO 

rn  100  IO*z |, NANGLE 

JCs20.*(?v-| . )+4 ,* (REG-2 . ) 

JCzJC*(ICA-l )»4 

CALL  PTA°F(  1  ,,’0,0,  10?4,,I0(  1  )) 

Fz(toA-n-icei?+i 

20  203  TzIFlL, IE1U 

X?z  F+I- 1  El L* 1 


I 
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1, 


I 


r 


i: 


ic 


c 


t 


JD(VD5  =  .!0(  I 5 
LA  i  *5  •*>  I CC5 
203  .'2(  LA  5 :  J0( !  +<005 

r 

C  CCNVFI,T  TO  CENTER  OF  »AFS  FYFTE" 

C 

.'0<Y5  =  JC<  55 
TWFTAiJOC 5) /ICO 
00  ICO  !=t,ICC5 
XF=X+I 

E.»F3rl04T(  I+ICCsj^FE 
IFCJOO'FI .NE.O)  00  TO  50 
JD(VF5r  10 
LB  =  VF+IC05 
JT(LP 5  =  1 00000 
50  COLL  FIMFl'P 

FCF: VF+I CC5 
vtu:  vc*;rf'5*icC5 
JC<XF5=  FLO<T(  .,2(*F5  5*p 

JE(wre):cLOAT(1tO<ve)  ) .  FLOAT!  J2(  YTF5  5/100. 

jdc/tn) =tc*  1  oo . 

100  CCNTIWF 

C 

C  OUTPUT  ON  T*°F:  *-FFCT10N  CN,Fo°0<:'*F^ ,  tnfta  or. 

c 

UPrPO.+CFUO- 1 .  5  +  < .  * ( °FG0-2 • 5 
.'??  =  YANGLE*2-I 
CALL  rT*pF( 2,UB,  1  ,  ITCH  ,U2(  1 5 5 
C 

C  CPE* TF  A°pAvc  *nd  DO  LEAPT  FOI'apFF  FIT 

C 

00  200  I  : 1,3072 
200  JO! 15=0 

20  30  L= ! El L, !E1U 
20  20  1=1, VANCLE 
X=L+! I  *  1 )*  ICC  1 2 
LC  :  y  -  !C26 
TMETO: FLC»T! JC!LC5 5/100. 

ELF2!I,  I  )=O0F!THETC»OA25 
L2  :  K  -  ICC7 

FLF2! 1,25 :FLOAT! J2!L255/10. 

LE  =  v  -  ICCF 
ELF 2 ( 1,35  =  FLOAT !  JI5!  LE5 5 
ELF2( 1 ,2  5: FLS2< 1,2  5*  FLF2( 1,2  5 
tF(FL?r<I,?5.fo.O.G5  CO  TO  30 
IF(ELFrcl,3).F0.0.05  CO  TC  30 
20  CONTINUE 

CALL  CBFI TF ( ELFT, A, XI , NPOL, NP VCLE 5 
25  IF  (ArATA.F'>.0.5  CO  TO  22 

C 

C  TTV  OUTPUT  CF  COF  ! 5 ,  FP°CP ,  y-SFCTI  0" ,  C»LC  POLY  X-FFOTI'H'. 
C 

IF(L  .E".  1 2 L"" 5  GO  TC  205 
GO  TO  22 

205  CA'.L  ITxpffL) 

IDU-  :  I  nr  +  10 
CALL  AFCIK212A.05 


/ 

£  , 


ro  21  10=1, mangle 

duelf3(io,  i> 

V2;  ELfr(I0,2> 
r3:ELSC( JC,3) 
n^ELrc<!o,4) 

CALL  TMY0UT(4,n,C2,D, 


RE4TE  OUTPUT  ARPAVS  (XI-XI,  TOTAL  X-SFCTICM,  Y-SFCTIO*'  F^C®) 


A4Pt4.*3.|415°*A( 1,1) 

LF  s  L  -  JFK.  *  1 
J0(LF):XI*1GGG. 

LG  :  L  -  rCC7  -I 
J0(LG):A4P/IG. 

LH  :  L  -  ICC6  -  1 
JO(LH):4< 1,2)*4G.*3.14I50 
DC  3G  Jts?,HPOL 
LlsL+ICC«-K.U-2)*ICC5 
L2iL*irn G-K  ,M  -2  >»  !OC5 
JCCL1  )i4(v»l,i  )» 10. *50000. 

J0(L2): 4(j 1 ,2)* 1G. 

CONTINUE 

C4LrUL4T!"r  OF  THE  RLCC*tfV  *  i>FC  \T.  FOP  OUTPUT  OF  arrays 
J?4  I  I  or  11/206  *1 
a3  :  UE  ♦  ,'B4 
4J?  i  J3 

ISV5  i  1.*  (»jb  *  4,)/2G. 

?W5  i  :?V5 

IR03  :  2.  *  ( AJP  -  20.  *  (SV5-l.))/4. 

RG5  r  IPGS 

JC5  s  20.  •  (CVS  -  1.)  *  4.  *  (P05  -  2.)  ♦  4. 

CALL  CT*PE(  2,  .'P5, 1 , 3072, U0(  1 ) ) 

OUTPUT  SV#  I  PEG4  POP  START  OF  OUTPUT  ARRAYS  -  3  PFCICNS  •  . 
CALL  • SO! !(?3, 5) 

A  JB3  r  UPS 

IPV6  t  l.*( aj35*4.)/20. 

SV6  i  IPV6 

IPG6  t  2 •  * ( » JP3-20 .* (SV6- I  .) )/4 . 

CALL  r  TVPE( I  SVG ) 

o*ll  ascii  (?s, 2) 

CALL  !TY°E< IPG6) 

CALL  ‘SC!!  (SI 0,2) 

LG1  :  IE1L-TC07-2 
CALL  IT^PEILGI) 

CALL  ASCII(P2G,1) 

00  TO  40 

erp 
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FUBPOUTIME  P3riT*(ELFC,A,*l  ,V?OL,VA'»CLn 
c  r»ra»4'«  accepts  a  vapiaflf  vpol  .if.*  a  pawl?  .lf.i?  . 

C  HPCL  IS  A  or  POL^OPIALS,  IT  1‘TPFASJVG  IT  CVF“  6,  VUFT  IVP«?AeF 

r  ri»”r  srLov.  vamglf  is  •  or  amglfs  or  imput  fata.  <nr  *c°  i5>. 
c  pfoc«am  "roiriEc  nior  csfit?  auc  7f,  ®FMAr?r  tve  fate. 

c  rtrcvrio'i  or  >«(>  fust  :  x**  wrrr 

C  PI VENE IOH  or  ELFC  4  A  KVFT  3E  FAFF  AF  IK  TVF  "» IV  P»OCPAT. 

c 

couplf  pwieion  H,L',rET,H?(j6) ,r?rrT 

DIPEKFICM  APG(ll>,rG(S),L(6),M(6),l'(l  I) 
rirtvsicv  aFC(2<i,A),A<6,2>,v(s,6) ,ppL'*m,e> 

F2UI  VALEVCE  <H2(  I  )*M<  I,  I )  > 

DAT*  FI  VCD, FI HCU2 /3M  SING,3MULAP  / 

vNCPKiNA'iGLE-KPOL 

CC  tt»  Jr  t ,  HA  V GIF 

apg(J)  i  n«r<j,n 

fAU  LEGPOl(FG,A»G( JJ.YPCL) 

DO  10  PsI.MWL 
10  POlYtJ,*)  I  Efi( *) 

c 

C  ®PLY  *««  FTCPEF  VALUE?  OF  PO  TM^OUCW  *<*JP('L-1>  AT  EAPV  OF  1?  avPLEF 
C  PCLYC  I  5,  VPCL ) 

C 


DC  20  Jr  I  ,'lPOL 
TO  20  *rl,MML 
?0  P<J,*>  i  0. 

PO  21  IrJ.'lPPL 
ro  21  Jr  I , VPOL 
rc  21  *r  1 , VA vQLt 

pi  wn,j>  r  ♦  pcLrtv.n.ppLro'.Ji/rLFrf 

CO  22  Jsl.MPOL 
PO  22  Irl.'lPCL 
K  r  1 ♦( J-l )*VP0L 
22  HP(V)  r  H(  I ,  J  ) 

PAIL  IMVV7(H2,s,P0L,rrT,L,'') 
irtCET)  27,25,27 
25  TAIL  ASCI  1 (f I KGU,2> 

RO  TO  100 


0 

P  V  IF  »»V  TWE*TFr,  CALCULATE  U 

c 

27  r  r 

DO  2J  J2  r  1,'JPOL 
ro  23  12  r  I.NPOL 
»  :  v-1 
I  r  VP0L-I2a1 
J  r  VPPL-J2H 
23  M<I,J>  r  M?(v) 

PO  30  Isl.HAVGLt 

30  I’CI)  r  0.0 

TO  31  1 r | , MPOL 
TP  31  Jr  I  ,  "A  'IGIE 

31  U<!)  :  l»(I)  ♦  ,*OL»(J,I  >»FtPr<J,J>/flFr< J,2) 


P 

P  PALPITATE  TVF  CPEFT*,  »I , 
C 


AND  P»Cff  FFPTICV 


»I?<'  s  0. 

re  j?  tn.'ipci 

J9  *<t, I)  :  0.0 

TO  40  Irl.VPOL 
ro  40  ,'r  I  ,  MPPL 

40  4<1,|>  :  4(1, I)  ♦  U(J)*H(J,I) 

C  NOV  ▼! 

C 

ro  45  iri,v#NGLE 
Jifiar  i  o. 
ro  44  ji i , v°ol 

44  fioar  :  FlGCt.0  +  4<J,  |>*?0LY<1,.I> 

EL?rn,4)  :  Pisrir 
ir<I.PT.M4!)PLE)  GO  TO  45 

ri  :  ELPra.j)  -  ELsr<i,4> 

»I  :  XI  •  VI 

vi ro  5  vi ro  ♦  »i/Ei?rc:t2) 

45  CONTINUE 

vi  :  vrpp/*»<(»pw 
TO  50  .'=I,"®01 

IT<H<v!,.t)  .LT.O.O)  M(J,J>  r  0.0 
4  ( J , .  )  t  r?ooT<H<.',J)> 

50  CO»'T!''UE 

c  *®og°4n  *rnic»-?  vi tm  cclculpte’'  vctct  in  Surrey),  rcrrrr 
r  4vr  £®®p°f  in  mi, .»),  4*ir  xi 

100  PETUP'i 

F‘»r 


c 

c 


<JO 


PUBPOl'TIVF  LECPOL  PECUPFIVE  EV4LU4TI0N  OF  LECFvrPF 

PO\.YNCTI4LF  ''POL  r  *  OF  POL YMOXI 4LF  TO  PETl'^N 
FUB90UT I NF  LEQPOL  <  "OL ,  »P0,  S'POL ) 

MFEMFION  POL(  6) 

XI  s  4 °G 
POL ( I )  s  1.0 
°0L<2)  s  XI 
ft  s  NPOL  -  I 
PO  09  U2,N 
FL  t  L 

POL(L*n  S  {<2.«F1-|  ,)*V|*PCL<L)  •  (*L-I  .)«?Ol(l-l  ))  /  ro 
PETUPN 

nr 
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C  5fB*0UTI"E  PINV7  —  T#*P!  FRO*  «L??ri'T!»,E  ri*iv  pF  jrt  «cir*T!nr 

C  ?UB*WT!'»S  P»C*»GE 

c  vote  ru-Pir  FPrr!ri("»  vbmbuff 

C  rr?C*IFT!r*!  «F  PBPB»ETEPS 

r  b  -  input  ■'PTRi*,  rrPT^OYEr  in  computation  bvo  ®f»<_»cfr 
c  Bv  fefultant  inverse 

c  n  -  c*nr®  of  pateif  * 

C  r  -  ’EFt’LTB'iT  Pf TEPFINBNT 

C  l  -  VCR*  VECTOR  OF  LENGTH  V 

r  X  -  VCPN  VECTOR  CF  LENGTH  “ 

c 

FL’FRCUTtVF  !MNV?<*,V,r,t,v) 

COUPLF  P*ECIfI0*'  », RIGA, woLP,C, PBBS 
DI»ENFIOH  B<  |),L(6),X<6) 
c  $SBPCU  pO®  LBPGE5T  ELEMENT 

0  s  1 .0000 
NX  I  -H 

ro  so  psi.n 

MV  s  -tr  ♦  M 
l<*>  i  * 

Pt«n  s  * 

RK  •  w  ♦  * 

RIGA  s  *<*v> 

DO  20  Ji*.  V 

r»  s  w*<j-n 

00  20  It*,M 
I J  *  !*  ♦  ! 

10  I F(  DBR?(?IGB>  -  Pa*>«<BCIJ>>  )  15,20,20 

15  PI  OB  r  »<!J> 

L<v>  :  I 
••(V)  r  J 
20  CONTINUE 

C 

C  INTERCHANGE  »0V« 

C 

J  •  L<F> 

!F<J-p>  55,25,25 
25  FI  :  *  •  1 

ro  so  t  s i , m 

PI  s  XI  ♦  V 
HOLP  :  -*(Pt) 

JI  s  »I  •  F  ♦  J 
A<*I>  r  A<JI> 

30  A<,« I>  *  HCLO 

C 


C  INTERCHANGE  COLUFNF 

0 

35  I  s  *<F> 

IF( !•*)  i5,<5,2B 
!•  JP  :  'P(I-I) 

PC  40  Jsl.N 
J*  I  W  ♦  J 
JI  s  JP  ♦  J 
“CLP  s 

*(JP)  s  A(Jt> 

bud  t  holt 


40 


PI  VIDE  CPU"'*!  5Y  NINUS  PIVOT  <  V*U‘E  CE  PIVOT  ELEMENT  t* 
CONTAINED  IV  PIGA  ) 

ir<pi(?*> 

D  :  O.OPOO 
petupm 
PO  95  Itl.N 
irn-Y>  50,55,50 
!»:»♦! 

AUK)  :  AUK)  /  (-BICA) 

CONTINUE 

nEtucr  katrix 

ro  69  I :  I ,  ‘I 
I*  :  NK  ♦  I 
HOLD  :  AUK) 

IJ  :  I  -  N 
DO  65  J*  I, S' 

IJ  J  IJ  *  1 

I'd-*)  60,65,60 

If<J-K>  62,55,62 

N  *  IJ  •  I  ♦  K 

AU»'>  i  MU.P *M*J)  *  HU) 

CONTINUE 

Pivicr  »0V  BY  PIVOT 

<rj  :  v  -  v 
TO  75  J:|,N 
KJ  :  «\l  ♦  N 
ir<J-Y>  70,75,70 
A<K4')  s  At* J)  /  913* 

CONTINUE 

PRODUCT  OE  •! VOTS 
p  i  r«sio* 

REPLACE  PIVfT  BY  RECIPRCm 

auk)  s  (i.oroo)  /  bio* 
continue 

EINAL  RCV  AND  CCIUPN  INTEPCHANCE 
*  s  M 

K  S  v  -  I 

IE(*)  190,190,105 

I  :  LC*> 

I  Ft  1**5  ISO, IPO, 10? 

,»C  s  '*<*-!> 

J*  s  M»U-I> 

PO  110  it:  I ,  V 
.**  :  J-  *  J 
NCIP  :  A(,'") 


JI  S  ,*•  ♦  J 

s  -MJl ) 

110  *(JI)  :  MOLD 

120  J  s  r<o 

If(J-K)  100,100,125 
125  *1  s  *  -  M 

P0  130  Is  1 , v 
ri  :  *!  *  M 
MOL 2  s  «<M!) 

JI  :  »l  •  »  ♦  J 

P<X!>  s  -»<v'I) 

130  *(JI)  r  MOLD 

00  TO  100 
150  METt’PN’ 

END 


oooooooo 
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APPENDIX  F 


Program  SIBPRP  and  Subroutine  HISTO 


cj  jsppo  e»f 

ITPIFirb  APRIL  l«7®  C»CN  \«LP°0,  ""cf.  ("IZ* LEY  4- 75 
CALCULATE?  ?P"nN  WlLl'iC  IV  PU"F  SILICON  TARGET? 
USING  PlICPVIllC  "FACTION, 

OPCTNATE  IS  CONCENT" ATI ON  f#?0»>C*V<A  OOPON  *  #SI)1  *FE 
A5SCISFA  I?  IN  "ICPONS  <I*X> 


CI^FV? !0V  NTSC  10241,  NTFArdOCA),  IPP"C( 102*) ,PPS<51 ."PTC?) 
corrov  NT 3,  NTS *v,  IPPPC 

C»TA  A|,  A2,A3,  A4/SHSIPP»,»y:»  ®CP,9HN)C  ,9W  4/7"/ 

FATA  P  1 , 32, ?3, P4/9H  IVPL',9WT  SAP.9MPLE  E.9HN  *  :/ 

FATA  Cl ,  02/9HPGN#  : ,  9M  rx:  / 

TATA  Dl,r?,DJ,r4/9H  1NPU.9MT  P0",5HCN  ,9«fY«:  / 

DATA  EI,G?,FJ/9M  0UT°,9MUT  F.9HV  •  :/ 

DATA  CCl,cr2,CC.J,rCA/0.00703«,  l?0.r?5e,33«.74«P,l  l«.?*9P/ 

I  CAUL  ASC!!<M,41 

c  input  Tirr  *r  "light  spectpup  prop  rri 
call  a«CTI(P|,4) 

CAUL  ANvT<{*<evi) 

CALL  AfCTICCI.il 
CAUL  ANYfNpC* Cn) 

C  INPUT  SA-PLE  OEPT“  INCRETENT  IN  PIC»ONS  CTT). 

CALL  ASCIKC2.I) 

CALL  ANYl^PCPr) 
r^rs/iocoo. 

,IS:?C.-<SWI.|  .)*4..(NGN-?.> 

C  INPUT  STANr»»C  (*U*E1  PORON  Tor  SPECTnU''  FTC*  PFCT»PF  l 
CALL  A  PC 11(01, 4) 

CAUL  ANvtNPCSWj) 
caul  Affiicci.n 

CALL  A»*vinP(»C*I) 

JLi20.»<fN!-l ,)»4.*(PCN-2.) 
c  INITIALISE  INPUT  an:  OUTPUT  ar»»v«  t?  ’F»C 
TO  too  1:1,1024 
NTSC I):0 
NTPA"<t):0 
I3P»0(I):0 
l 00  CONTINUE 

CAUL  OTAAFC  l,.'f,0,l024.NT'A««<  m 
CALL  CTAPFC  I,.'L,0,I024,NTFC  ID 

C  INPUT  TI-E  cr  "LIGHT  "PECT’A  SHOULD  PE  LAPFLFD  VITV  RUNNING 
C  PAPANFTr»£  Af  POUUOVSt  CM  0:CHANNrU  VIOTM  TM  *?EC:  C«  l  : 
c  QACPA  PA*  CHANNEL  TT"ES  TEN:  CH  ?:  FLICVT  ®ATN  In  Cr:  C«  *: 

C  INTEGRATED  PPCTPN  PEA*  CURRENT  IN  ICPUL:  CH A: P"C TON  E?.’E°PY  "EV  *  lt-f.0 
ro  200  1:1,9 

•ppci>:'vrrc!> 

RPf  (  I ):*:TS»P< I ) 

200  ISPPOC I):NTf A»( I) 


calculatf  tf*o  nrr  for  eac*  tc"  rprcTFL*1 

?P:*P?C2)/I0.*JJ.4«PP?C?)/«PS(I) 
•*EA'':®«EC2>/tt.*Ji.4«RPF(J)/PPS(  11 
FP":«P«( 91/1000. 

srr  ccmc entpatic**  e^ual  tc.  -rpc. 

•:0. 

SAHPUt  NU*?ER  OENSIT"  "OR  «ILICCM  <?0»n*! 


12?. I*  E*T 1 1 


*y 


88 


301 

C 


190 

C 

c 

c 

203 

C 

c 

c 

220 

209 

20* 

209 

C 

206 

f 


C 

2i0 

C 


c 

2*0 

c 


208 

C 

20® 

300 

r 

301 


SN0:A*.®7E*2l 

START  ENERGY  OF  ?®OTPN  TEN  STEPS  CFRONT  C*  TA»CFT. 

PO  190  T:  I  ,  10 

EPSAP*J.4Plf-2l*<AL0G(  EPM  )■*? .3 1«>/FPH 

f®h*  eph+snt*  c  eps  a>*>  -rx 

CALCULATE  CONCENT® ATI Cm  *S  A  FUNCTION  CF  2FPT«  (T») 

PC  300  Is  II, 1024 

CALCULATE  STO ®PING  PONE?  EO®  ?C?C*’  *  SILICON 

“JOT  °EP  IIP  ATCY  *S  WAS  CCNE  IN  r-SFOT. 

EPSP: I . 324 E-2 I *< AL0G(EPH)*2.® 69 )/EP« 

EPSA*'i2  ,4?1E*2 1  *<  *LOG(  EPH)*2.3  1®)/*?H 
EPS«A:EPSB/EPS*r 

CALCULATE  TVS  NFV  rNEPG*  OF  TUE  P*CTPV  ®ASE2  ON  THE 
CONCENTRATION  OF  THE  BORON  PRESENT  USING  CURRENT  CONCENTRATION 
ANP  STOPPING  ROWER. 

IFCR-I.)  203,204,204 
EP’.:FPH-SNC«(EPSAr*R*EPS?/<l  .O-®))*!"* 


GO  TO  209 

EPL*  EPM- 128. 14E+21 *EPS®*ry 
IF(EPL-3.0P)  301 ,206,206 

CALCULATE  NEUTRON  ENE®Gw  r«0T  ''I'TTATICS  EC?  0  TEC  i  PI VF*  r® 
ENL:CCl»EPL*<CC2-rC3/EPL*'2.*S^»T(CC4-rC3/EPL>) 

EN*.'i  CO  I  ■EPH*(CC2*CC3/EPH+'2.»SO?T(CC4-CC9/rPR) ) 

CALCULATE  NEUTRON  OUANNEL  *  <*FLI TI VISTICALL*)  EO®  S»“PLF 
TLSsO,  7235*  »PS< 3 )  /SC pT<  ENL  ) 

TUS:0.7236«»PS<3)/SPRT(ENL') 

PlLSr’RA*-TlE*IOOO./*PSC  I ) 

C  IUSr  ’ SA* - T US*  1 000 . /RPS (  1  ) 

I r(C ILS)  301,301,230 

CALCULATE  *  OE  COUNTS  »ETWEEN  GIVEN  CHANNELS,  EO*  SAMPLE  SPECT®l,v 
CALL  “! <“T0<  OILS «  OIL'S, NTS AP,  SIP'S*  v) 

CALCULATE  NEUTRON  C« anuEL  »  <*ELI  TIVISTIOALL'') ,  EC?  ?0e0N  SPE CTRL'S 
TL5: 0. 7236*®P9<  3) /  S^RT (  ENL ) 

TU?i0.7236*?P8<3)/SCRT<FNU> 

CIL?=  ?3*TLS* 1000 , /PPP( I ) 
riUBrTE-TU?»IOOO./RP9< 1) 

IE  tCILT)  301,301.2*0 

CALCULATE  •  OE  COUNTS  ®ETVE EN  GIVEN  CHANNELS,  FOR  ?0»0N  ?PFPTBIT 
CALL  “ISTP<  CILE,  Crt.'?,NTP,  SUSP  ) 

CALCULATE  CONCENTPATION 

?0r 0P2( A)»  FPSA v* SUrrAr/(?»«<3)*RPR(3) ) 
R2iRPS(4)*EPS?*SUTP/(?PS(3)*RPS(3)> 

P I  r  ®P/ ( *o« PC* FPfRA  *?2) 

IE  (PI)  1,20®, 20« 

If(E I .CT.1.0)  *1*1,0 
rc^Ps (®i-p)/?i 


RrPt 

CHECF  TO  PE  Fl'RE  THAT  NEW  R  IS  VIT«IN  I*  CE  T«F  OLT  " 
I E( ( CPFP ,GT ,0,01 ) .OR .  ( OC'P.LT .*0 .0 1 ) )  00  TO  220 
IPP*0<  f)*R* I ,OE*S 
EPHiEPL 
CALL  LITES(I) 

CCKT I  NOT 


OUTPUT  SW  4  RON  ON  CTI 
CALL  ASCIt(rt,3) 
CALL  ANyInp<?v!) 
CALL  ASri!cri,l) 


c>  a’ 


Nvv 


•r*  'nJ 


call 

JCi20.«<  <^1-1  .5*4  .  •(  ?GN'-2  . ) 

call  :t4pe<i,jc,i,io?4,i?pr(?(|)) 

OUTPUT  STAPTS  IN  SPA  CHANNEL  10.  FFCNT  SURFACE  0r  TA°GrT  Tc 
N0HINALLv  IN  PPA  CHANNEL  20.  0PCI',ATE  IS  CCHC*‘’TPATIo\r  n*re  r*» 
GO  TC  1  '  " 

STOP 
ENO 


SL’2 eOl’T I "F  WI STO( CIL 
?1 TENSION  *TOF(  1 024) 

11  "  CIL*. 5 

12  =  CIL*. 5 
CELTA  i  Cl li- CIL 
P  :  I?-II 

sur  i  0.0 
I  E(  1 1 )  1  Ot  1 0, 1 
ir<TELTA-I .52,2,3 

I  ECO  4, 4 ,3 

I“C  T-  1 54,  7, 5 
PI  =  CELTA 
re  r  0.0 
GO  TO  ? 

VI  s  11*1 

re  i  12- i 

20  S  LsNI,'2 

ICV  r  L* I 

*C»TA  :  'ITOF(ICH) 

sur  i  sur* re a ta 

*11  s  II 

*12  i  12 

FI  r  0.5-(CIL-*It> 

F2  t  0.5*01 ’J-TI2 
irw  S  11*1 
VSATA  r  NTOF( ICu) 

SU"  r  sur  *  FI  ■  VSATA 

ICV  -  J2*| 
re*T4  r  VTCOfTrn) 

sur  :  su“*F2*rrAT« 


90 


REFERENCES 

1.  F.  Ajzenberg-Selvove ,  Nuclear  Physics  A  248.  i  (1975). 

2.  J.  P.  Blaser,  F.  Boehm,  P.  Marmier  and  P.  Scherrer,  Helv.  Phvs. 

Acta  24,  465  (1951). 

3.  J.  K.  Bair,  J.  D.  Kington  and  H.  B.  Willard,  Phys.  Rev.  100. 

21  (1955). 

4.  J.  H.  Gibbons  and  R.  L.  Macklin,  Phys.  Rev.  114,  571  (1959), 

5.  J.  B.  Marion,  T.  W.  Bonner  and  C.  F.  Cook,  Phys.  Rev.  100.  91  (1955). 

6.  G.  J.  F.  Legge  and  I.  F.  Bubb,  Nucl.  Phys.  26^  616  (1961). 

7.  R.  E.  Segel,  S.  S.  Hanna  and  R.  G.  Allas,  Phys.  Rev.  139.  B818 
(1965). 

8.  J.  C.  Overley  and  R.  R.  Borchers,  Nucl.  Phys.  65,  156  (1965). 

9.  L.  Van  der  Zwan  and  K.  W.  Geiger,  Nucl.  Phys.  A  306,  45  (1978). 

10.  J.  C.  Overley,  R.  P.  Ebright,  and  H.  W.  Lefevre,  IEEE  Transac¬ 
tions  on  Nuclear  Science  NS-26.  1624  (1979). 

11.  C.  A.  Burke,  Ph.  D.  Thesis,  University  of  Oregon,  1972  (unpub¬ 
lished)  . 

12.  J.  C.  Overley,  private  communication. 


91 


13.  J.  F.  Janni,  Air  Force  Weapons  Laboratory  Tech.  Report  No. 
AFWL-TR-150  (1966)  (unpublished). 

14.  E.  Segre,  Experimental  Nuclear  Physics,  John  Wiley  4  Sons,  Inc., 
New  York,  Vol.  1,  168  (1953). 

15.  Weast,  Handbook  of  Chemistry  and  Physics,  52nd  Ed.,  Chemical 
Rubber  Co.,  Ohio,  F  18  (1972). 

16.  M.  T.  Lunnon,  Ph.  D.  Thesis,  University  of  Oregon,  1974  (unpub¬ 
lished)  . 

17.  P.  R.  Bevington,  W.  W.  Roland,  and  H.  W.  Lewis,  Phys.  Rev.  121, 
871  (1961). 

18.  R.  R.  Borchers  and  C.  H.  Poppe,  Phys.  Rev.  129,  2679  (1963). 


19.  H.  Liskien  and  A.  Paulsen,  Atomic  Data  and  Nuclear  Data  Tables, 
15-1,  58  (1975). 


20.  H.  W.  Lefevre,  The  Absolute  Neutron  Spectrum  to  E  *  2.4  MeV 


from  the  (p,n)  Reaction  in  Thick  Lithium  Metal  Target,  (unpublished). 


21.  D.  J.  Hughes  and  R.  B.  Schwartz,  Neutron  Cross  Sections,  2nd  Ed., 
Associated  Universities,  Inc.,  Brookhaven  National  Lab,  New  York, 
38  (1958). 


22.  J.  C.  Overley  and  H.  W.  Lefevre,  Proc.  Symposium  Radiation  Effects 
on  Solid  Surfaces.  N.  Kaminsky,  Ed.,  American  Chemistry  Society, 
Washington,  D.  C.,  282  (1976). 


92 

23.  J.  C.  Overley  and  H.  W.  Lefevre,  Scientific  and  Industrial 
Applications  of  Small  Accelerators.  Fourth  Conference-1976, 

J.  Duggan  and  I.  Morgan,  Ed.,  IEEE,  New  York,  NY,  76  CH  1175-9NPS, 
557  (1976). 

24.  C.  A.  Bams,  J.  C.  Overley,  2.  E.  Switkowski,  and  T.  A.  Tombello, 
App.  Phys.  Lett.,  31,  239  (1977). 

25.  C.  W.  White,  J.  Naravan,  R.  T.  Young,  Science  204 ,  461  (1979). 

26.  K.  Wittmaack,  J.  Maul  and  F.  Schulz,  Ion  Implantation  in  Semi¬ 
conductors  and  Other  Materials,  B.  L.  Crowder,  Ed.,  Plenum  Press, 
New  York,  147  (1973). 

27.  D.  Lecrosnier,  J.  Paugam  and  J.  Gallou,  Appl.  Phys.  Letters,  30, 
323  (1977). 

28.  D.  Eirug  Davis,  Appl.  Phys.  Letters,  .13,  243,  (1968). 

29.  D.  Eirug  Davis,  Can.  J.  Phvs.,  47_,  1750  (1969). 

30.  H.  Ryssel,  H.  Kranz,  K.  Muller,  R.  Henkelmann  and  J.  Biersach, 
Appl.  Phys.  Letters,  30,  399  (1977). 

31.  Y.  Ohmura,  K.  Koike  and  H.  Kobayashi,  Ion  Implantation  in  Semi¬ 
conductors  ,  Susumu  Namba,  Ed.,  Plenum  Press,  New  York,  183  (1974). 


32.  Y.  Akasaka  and  K.  Horie,  Ion  Implantation  in  Semiconductors  and 
Other  Materials,  B.  L.  Crowder,  Ed.,  Plenum  Press,  New  York, 

183  (1974). 


£ 

& 


